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(N ' Abstract: This course provides a self contained introduction to the general 

theory of rclativistic brane models, of the category that includes point particle, 
string, and membrane representations for phenomena that can be considered as 
£Ni ■ being confined to a worldsheet of the corresponding dimension (respectively one, 

t**"* \ two, and three) in a thin limit approximation. The first part of the course is 

concerned with purely kinematic aspects: it is shown how, to second differential 
order, the geometry (and in particular the inner and outer curvature) of a brane 
worldsheet of arbitrary dimension is describable in terms of the first, second, 
| and third fundamental tensor; the extension to a foliation by a congruence of 

such worldshccts is also briefly discussed. In the next part, it is shown how 
- to lowest order in the thin limit - the evolution of such a brane worldsheet 
will always be governed by a simple tensorial equation of motion whose left 
(\y\ hand side is the contraction of the relevant surface stress tensor T^ v with the 

i (geometrically defined) second fundamental tensor K^ u p , while the right hand 

side will simply vanish in the case of free motion and will otherwise be just the 
orthogonal projection of any external force density that may happen to act on 
r% \ the brane. (Allowance for first order deviations from such a thin limit treat- 

jj^ ■ ment would require evolution equations of a more complicated kind of which 

a prototype example is presented.) The last part of the course concentrates 
on the case of a string: in a four dimensional spacetime background, this is 
only case that is non trivial in the strong sense of having a worldsheet with 
both dimension and codimension greater than one. This case is of particular 
importance for cosmology because it applies to the cosmic strings that arise 
as vortex defects of the vacuum in many kinds of field theory, and that may 
conceivably play an important role in the evolution of the universe. It is shown 
how to set up a Lagrangian providing the complete system of (both internal and 
external) string evolution equations governing the large scale motion of vacuum 
vortices of both the simple Nielsen-Olesen-Kibble type and of Witten's super- 
conducting type, including allowance for conservative external forces of both 
electromagnetic and Kalb-Ramond axion type, but neglecting the dissipative 
effects that are likely to be important soon after string formation but not later 
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on. The last part of the course deals with closed loops under typical circum- 
stances in which not only dissipative but also conservative (electromagnetic or 
axion type) external forces can be neglected in a lowest order approximation, 
so that the evolution equations will be homogeneous, and thus easily soluble 
both for dynamic circular configurations and also for the stationary (centrifu- 
gally supported) configurations known as vortons, which, if they arc sufficiently 
stable, may be of considerable cosmological significance. 



As an aide memoire for readers equiped for viewing colour, warmer hues will 
be used for physical quantities of a material nature, as contrasted with essen- 
tially geometric quantities for which cooler hues will be reserved. Among the 
latter scalar and tensorial fields determined just by the background spacetime 
geometry itself will be indivated by blue symbols while background coordinate or 
frame dependent quantities will be indicated by black symbols. Similarly scalar 
and tensorial fields determined on the brane worldsheet just by the embedding 
geometry will be indicated by green symbols quantities that are fully determined 
as scalars or tensors just by the embedding geometry, while internal coordinate 
or frame dependent quantities on the brane worldsheet will be indicated by sky- 
blue symbols. As far as more essentially physical quantities are concerned, pink 
symbols will be used for background gauge fields while magenta symbols will be 
used for physical background fields. On the other hand beige symbols will be 
used for worldsheet supported gauge fields while brown symbols will be used for 
physical quantities confined to and intrinsic with respect to the worldsheet, and 
violet symbols will be used for quantities of mixed type, including physical fields 
confined to the worldsheet but defined with respect to an external reference sys- 
tem, as well as external fields that may be defined with respect to fields on the 
world sheet. Finally red symbols will be used for global integrals and coupling 
constants. 
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0.1 Introduction 



The present course is an updated and extended version of lectures written up 
a couple of years ago for a summer school on "Formation and Interactions of 
Topological Defects" jjj at the Newton Institute in Cambridge. In preparation 
for the more specific study of strings in the later sections, the first part of this 
course is intended as an introduction to the systematic study, in a classical rela- 
tivistic framework, of "branes" , meaning physical models in which the relevant 
fields are confined to supporting worldsheets of lower dimension than the back- 
ground spacetime. While not entirely new ||, ||, this subject is still at a rather 
early stage of development (compared with the corresponding quantum theory 
which has been stimulated by the rise of "superstring theory"), the main 
motivation for recent work Q on classical relativistic brane theory being its ap- 
plication to vacuum defects produced by the Kibble mechanism j(| , particularly 
when of composite type as in the case of cosmic strings attached to external 
domain walls [Q and of cosmic strings carrying internal currents of the kind 
whose likely existence was first proposed by Witten ||. A propos of the latter, 
a noteworthy example of the progress that has been achieved since the prev ious 
version Q of these notes was written is the construction (see Subsection 5.7) for 
the first time of a class of reasonably simple analytically explicit string models 
H that arc capable of describing the macroscopic behaviour of such Witten 
vortices in a manner that is qualitatively satisfactory and quantitatively ade- 
quate not just when the current is very small but even for the largest physically 
admissablc values. This makes it possible to provide an analytic description 
(see Subsection 3.5) of the dynamics of the circular ring configurations whose 
equilibrium states (see Subsection 3.6) are the simplest examples of what are 
known as vortons. 

Before the presentation in Section | of the dynamic laws governing the evo- 
lution of a brane worldsheet, Section [l provides a recapitulation of the essential 
differential geometric machinery jt(^, |fl| needed for the analysis of a timelike 
worldsheet of dimension d say in a background space time manifold of dimension 
n. At this stage no restriction will be imposed on the curvature of the metric - 
which will as usual be represented with respect to local background coordinates 
x^ 1 (fi~ 0, n-1) by its components g^ v - though it will be postulated to be 
flat, or at least stationary or conformally flat, in many of the applications to be 
discussed later. 



1 Worldsheet curvature analysis 
1.1 The first fundamental tensor 

The development of geometrical intuition and of computationally efficient meth- 
ods for use in string and membrane theory has been hampered by a tradition of 
publishing results in untidy, highly gauge dependent, notation (one of the causes 
being the undue influence still exercised by Eisenhart's obsolete treatise "Rie- 
mannian Geometry" Il3]). For the intermediate steps in particular calculations 
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it is of course frequently useful and often indispensible to introduce specifically 
adapted auxiliary structures, such as curvilinear worldshcet coordinates a 1 (i= 
0, d-1) and the associated bitensorial derivatives 

or specially adapted orthonormal frame vectors, consisting of an internal subset 
of vectors t A M ( a= 0, ... , d-1) tangential to the worldsheet and an external 
subset of vectors X x fl ( x = 1, — , n-d) orthogonal to the worldsheet, as 
characterised by 

A^1a B , t/A xM = 0, K^Kn = S XY , (2) 

where r\ AB is a fixed d-dimensional Minkowski metric and the Kronecker matrix 
S XY is a fixed (n-d)-dimensional Cartesion metric. Even in the most recent 
literature there are still (under Eiscnhart's uninspiring influence) many exam- 
ples of insufficient effort to sort out the messy clutter of indices of different 
kinds (Greek or Latin, early or late, small or capital) that arise in this way by 
grouping the various contributions into simple tensorially covariant combina- 
tions. Another inconvenient feature of many publications is that results have 
been left in a form that depends on some particular gauge choice (such as the 
conformal gauge for internal string coordinates) which obscures the relationship 
with other results concerning the same system but in a different gauge. 

The strategy adopted here [jl3| aims at minimising such problems (they can 
never be entirely eliminated) by working as far as possible with a single kind 
of tensor index, which must of course be the one that is most fundamental, 
namely that of the background coordinates, x^. Thus, to avoid dependence on 
the internal frame index a (which is lowered and raised by contraction with 
the fixed d-dimcnsional Minkowski metric rj AB and its inverse rj AB ) and on the 
external frame index x (which is lowered and raised by contraction with the 
fixed (n-d)-dimcnsional Cartesian metric S XY and its inverse 8 XY ), the separate 
internal frame vectors t A M and external frame vectors A Y M will as far as possible 
be eliminated in favour of the frame gauge independent combinations 

< = ^^,, J* (3) 

of which the former, rj^, is what will be referred to as the (first) fundamental 
tensor of the metric, which acts as the (rank d) operator of tangential projection 
onto the world sheet, while the latter, A?„, is the complementary (rank n-d) 
operator of projection orthogonal to the world sheet. 

The same principle applies to the avoidance of unnecessary involvement of 
the internal coordinate indices which are lowered and raised by contraction with 
the induced metric on the worldsheet as given by 

7ij = g^x^iX^j , (4) 

and with its contravariant inverse 7 y . After being cast (by index raising if neces- 
sary) into its contravariant form, any internal coordinate tensor can be directly 
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projected onto a corresponding background tensor in the manner exemplified 
by the intrinsic metric itself, which gives 

rT = 7 y <,< J - > (5) 

thus providing an alternative (more direct) prescription for the fundamental 
tensor that was previously introduced via the use of the internal frame in (||). 
This approach also provides a direct prescription for the orthogonal projector 
that was introduced via the use of an external frame in (|^) but that is also 
obtainable immediately from (^J) as 

-^ = <?t-<- (6) 
As well as having the separate operator properties 

rf p rf v = rf vi £ p £ v = J? v (7) 

the tensors defined by (|J) and (^) will evidently be related by the conditions 

V ^ v = = £ p rf v . (8) 



1.2 The inner and outer curvature tensors 

In so far as we are concerned with tensor fields such as the frame vectors whose 
support is confined to the d-dimensional world sheet, the effect of Ricmannian 
covariant differcntation V M along an arbitrary directions on the background 
spacetime will not be well defined, only the corresponding tangcntially projected 
differentiation operation 

%=V v ^u, (9) 

being meaningful for them, as for instance in the case of a scalar field ip for 
which the tangcntially projected gradient is given in terms of internal coordinate 
differentiation simply by V"^ = j^x^^ ip t ji . 

An irreducible basis for the various possible covariant derivatives of the frame 
vectors consists of the internal rotation pseudo-tensor p " and the external 
rotation (or "twist" ) pseudo-tensor w ■ 1 ' as given by 

P»"p = rfo i A P % i A ° = -p» p v , = A, ^ p%K- ° = -V , (10) 

together with their mixed analogue K pv p which is obtainable in a pair of equiv- 
alent alternative forms given by 

K^p = J? CT i A a = -rf v \ x p %\ x « . (11) 

The reason for qualifying the fields (|l0|) as "pseudo" tensors is that although 
they are tensorial in the ordinary sense with respect to changes of the back- 
ground coordinates x^ they are not geometrically well defined just by the geom- 
etry of the world sheet but are gauge dependent in the sense of being functions 
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of the choice of the internal and external frames t A ^ and A Y M . The gauge de- 
pendence of p". and 'oo^p means that both of them can be set to zero at any 
chosen point on the worldsheet by choice of the relevant frames in its vicin- 
ity. However the condition for it to be possible to set these pseudo-tensors to 
zero throughout an open neigbourhood is the vanishing of the curvatures of 
the corresponding frame bundles as characterised with respect to the respec- 
tive invariance subgroups SO(l,d-l) and SO(n-d) into which the full Lorcntz 
invariance group SO(l,n-l) is broken by the specification of the d-dimensional 
world sheet orientation. The inner curvature that needs to vanish for it to be 
possible for p" to be set to zero in an open neighbourhood is of Riemannian 
type, is obtainable (by a calculation of the type originally developed by Cartan 
that was made familiar to physicists by Yang Mills theory) as Jl(| 

RkX^u = 27A ?? > W [A V K]P / r + 2p [K ^p x]nv , (12) 

while the outer curvature that needs to vanish for it to be possible for the "twist" 
tensor to be set to zero in an open neighbourhood is of a less familiar type 
that is given by 

= 2J? ff J_ T M rr [x V K] wS T + 2w [K ' m w x]m . (13) 



The frame gauge invariance of the expressions ( |12| ) and (13), - which means 
that Rkx^v and Okx^v are unambiguously well defined as tensors in the strictest 
sense of the word - is not immediately obvious from the foregoing formulae, but 



it is made manifest in the the alternative expressions given in Subsection 1.5 



1.3 The second fundamental tensor 

Another, even more fundamentally important, gauge invariance property that is 
not immediately obvious from the traditional approach - as recapitulated in the 
proceeding subsection is - that of the entity K ^ v p defined by the mixed analogue 
( pr] ) of ([To|), which (unlike p u and w ', but like RkX^v and i? K A M i/) is m fact a 
geometrically well defined tensor in the strict sense. To see that the formula ( pi] ) 
does indeed give a result that is frame gauge independent, it suffices to verify 
that it agrees with the alternative - manifestly gauge independent definition [j5| 

*V= ri a „%V p a - (14) 

whereby the entity that we refer to as the second fundamental tensor is con- 
structed directly from the the first fundamental tensor rf as given by (|5|). 

Since this second fundamental tensor, K Aljy p will play a very important role 
throughout the work that follows, it is worthwhile to linger over its essential 
properties. To start with it is to be noticed that a formula of the form (|l4|) could 
of course be meaningfully meaningful applied not only to the fundamental pro- 
jection tensor of a d-surface, but also to any (smooth) field of rank-d projection 
operators n ^ v as specified by a field of arbitrarily orientated d-surface elements. 
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What distinguishes the integrable case, i.e. that in which the elements mesh to- 
gether to form a well defined d-surface through the point under consideration, is 
the condition that the tensor defined by (|l4|) should also satisfy the Weingarten 
identity 

K^f = (15) 

(where the square brackets denote antisymmetrisation) , this symmetry prop- 
erty of the second fundamental tensor being derivable §, |oJ as a version of 
the well known Frobenius theorem. In addition to this non-trivial symmetry 
property, the second fundamental tensor is also obviously tangential on the first 
two indices and almost as obviously orthogonal on the last, i.e. 

r fl K a „P=K fl „°r )a p = 0. (16) 

The second fundamental tensor has the property of fully determining the 

tangential derivatives of the first fundamental tensor r\ ^ v by the formula 

= 2K tl ( yp ) (17) 

(using round brackets to denote symmctrisation) and it can be seen to be char- 
acterisable by the condition that the orthogonal projection of the acceleration 
of any tangential unit vector field will be given by 

= J^it" , ii^ = u v V u u» . (18) 



1.4 The extrinsic curvature vector and the conformation 
tensor 

It is very practical for a great many purposes to introduce the extrinsic curvature 
vector K^, defined as the trace of the second fundamental tensor, which is 
automatically orthogonal to the worldsheet, 

= K v v » , rf v IC = . (19) 

It is useful for many specific purposes to work this out in terms of the intrinsic 
metric 7^ and its determinant |7|. It suffices to use the simple expression 
W = y^x^^ipj for the tangentially projected gradient of a scalar field Lp on the 
worldsheet, but for a tensorial field (unless one is using Minkowski coordinates 
in a flat spacetime) there will also be contributions involving the background 
Riemann Christoffel connection 

= 9 va {g^, P ) - • (20) 

The curvature vector is thus obtained in explicit detail as 

K» = %rT = -if (vlThXi) „■ + 7 y <i^V p . (21) 



This last expression is technically useful for certain specific computational pur- 
poses, but it must be remarked that much of the literature on cosmic string 
dynamics has been made unnecessarily heavy to read by a tradition of working 
all the time with long strings of non tensorial terms such as those on the right of 
( f2lj ) rather than taking advantage of such more succinct tensorial expressions as 
the precccding formula V^?/^. As an alternative to the universally applicable 
tensorial approach advocated here, there is of course another more commonly 
used method of achieving succinctness in particular circumstances, which is to 
sacrifice gauge covariance by using specialised kinds of coordinate system. In 
particular for the case of a string, i.e. for a 2-dimensional worldshcct, it is stan- 
dard practise to use conformal coordinates a a and c 1 so that the corresponding 
tangent vectors ± M = x p a and x'^ = x^ satisfy the restrictions x p x' ^ = 0, 
x^Xfj, + x' p x' M — 0, which implies \/||7|| = x'^x' M = —x^x^ so that (|T]) simply 
gives y/f^IC = x" v -x v + {x'^x'p -i^)ry p . 

The physical specification of the extrinsic curvature vector (jllj) for a time- 
like d-surface in a dynamic theory provides what can be taken as the equations 
of extrinsic motion of the d-surface || , , the simplest possibility being the 
"harmonic" condition K p = that is obtained (as will be shown in the fol- 
lowing sections) from a surface measure variational principle such as that of 
the Dirac membrane model [^J, or of the Goto-Nambu string model || whose 
dynamic equations in a flat background arc therefore expressible with respect 
to a standard conformal gauge in the familiar form x"^ — 'x p = 0. 

There is a certain analogy between the Einstein vacuum equations, which 
impose the vanishing of the trace 1Z^ V of the background spacctime curvature 
TZ-\n P v, and the Dirac-Gotu-Nambu equations, which impose the vanishing of 
the trace K v of the second fundamental tensor K\^ u . Just as it is useful to sep- 
arate out the Weyl tensor |Q, i.e. the trace free part of the Ricci background 
curvature which is the only part that remains when the Einstein vacuum equa- 
tions are satisfied, so also analogously, it is useful to separate out the the trace 
free part of the second fundamental tensor, namely the extrinsic conformation 
tensor ]Tc| |, which is the only part that remains when equations of motion of 
the Dirac - Goto - Nambu type are satisfied. Explicitly, the trace free extrinsic 
conformation tensor C \ IV P of a d-dimensional imbedding is defined |ic| in terms 
of the corresponding first and second fundamental tensors rj^ and K ^ v p as 

C^" = K^p - ^n^K" , CV - . (22) 

Like the Weyl tensor W\^ p y of the background metric (whose definition is given 



implicitly by (27) below) this conformation tensor has the noteworthy property 
of being invariant with respect to conformal modifications of the background 
metric: 

(23) 

This formula is useful for calculations of the kind undertaken by Vilenkin 
fp8| in a standard Robertson- Walker type cosmological background, which can 
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be obtained from a flat auxiliary spacetime metric by a conformal transformation 
for which e Q is a time dependent Hubble expansion factor. 

1.5 *The Codazzi, Gauss, and Schouten identities 

As the higher order analogue of (|l4|) we can go on to introduce the third funda- 
mental tensor H as 

Zx^ p = V > \^xK aT a , (24) 

which by construction is obviously symmetric between the second and third 
indices and tangential on all the first three indices. In a spacetime background 
that is flat (or of constant curvature as is the case for the DeSittcr universe 
model) this third fundamental tensor is fully symmetric over all the first three 
indices by what is interpretable as the generalised Codazzi identity which is 
expressible JTo[ in a background with arbitrary Riemann curvature 'R.\ ll p a as 

Exp," = ^(a m ,) p + \if^\^ a v) Tl a / a ^ (25) 

It is to be noted that a script symbol 1Z is used here in order to distinguish 
the (n - dimensional) background Riemann curvature tensor from the intrinsic 
curvature tensor ([l^) of the (d - dimensional) worldship to which the ordinary 
symbol R has already allocated. 

For many of the applications that will follow it will be sufficient just to treat 
the background spacetime as flat, i.e. to take 7^ CTT ^ Q = 0. At this stage however, 
we shall allow for an unrestricted background curvature. For n> 2 this will be 
decomposible in terms of its trace free Weyl part W A1 ^ P (T (which as remarked 
above is conformally invariant) and the corresponding background Ricci tensor 
and its scalar trace, 

1l pv = ftp/, , K = K V V , (26) 

in the form Q 

(in which the Weyl contribution can be non zero only for n> 4). In terms of 
the tangential projection of this background curvature, one can evaluate the 
corresponding internal curvature tensor ( [l2| ) in the form 

R MV " a = 2K" ] /K V] „ + V > x „1l KX a T r, " aV \ , (28) 

which is the translation into the present scheme of what is well known in other 
schemes as the generalised Gauss identity. The much less well known analogue 
for the (identically trace free and conformally invariant) outer curvature ( p"3| ) 
(for which the most historically appropriate name might be argued to be that 
of Schouten |l3) is given Jlp[] in terms of the corresponding projection of the 
background Weyl tensor by the expression 

= 2C { ^ p C v]Ta + V > \W K X a rV a ■ (29) 
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It follows from this last identity that in a background that is flat or conformally 
flat (for which it is necessary, and for n> 4 sufficient, that the Weyl tensor 
should vanish) the vanishing of the extrinsic conformation tensor C fIl/ p will be 
sufficient (independently of the behaviour of the extrinsic curvature vector K^) 
for vanishing of the outer curvature tensor flp, v p which is the condition for 
it to be possible to construct fields of vectors A M orthogonal to the surface and 
such as to satisfy the generalised Fermi- Walker propagation condition to the 
effect that Jf^V^Ap should vanish. It can also be shown (l(J (taking special 
trouble for the case d=3 ) that in a conformally flat background (of arbitrary 
dimension n) the vanishing of the conformation tensor C^ p is always sufficient 
(though by no means necessary) for conformal flatness of the induced geometry 
in the imbedding. 

1.6 *Extension to a foliation 

It is useful for many purposes to extend the analysis of a single embedded d- 
surface (as described in the preceeding subsections) to the more general context 
of a smooth foliation by a congruence of such surfaces. Even when the object of 
physical interest is just a single embedded surface it is geometrically instructive 
to consider its relation to neighbouring surfaces of a purely mathematical nature, 
and of course there are also many applications involving modes in which a 
foliation by congruence of surface is a physically essential constituent structure 
(for example the congruence of flow trajectories in an ordinary fluid model ]l5| , 
or the congruence of string- like vortex lines in the kind of model [ fl6| appropriate 
for the macroscopic description of a rotating superfmid.) 

In the case of such a foliation, the first fundamental tensor rj^ as specified 
by (|5|) and of the corresponding orthogonal projector Jt v as specified by (^) will 
no longer have support confined to a single d-surface, but will be well defined 
as ordinary tensor fields over at least an open neighbourhood of the background 
spacctimc. This means that they will have well defined (Ricmannian or pscudo- 
Riemannian) covariant derivatives of the ordinary unrestricted kind (not just 
tangential covariant derivatives of the kind described in the preceeding sections). 
These covariant derivatives given will be fully determined by the specification 
of a certain (first) deformation tensor, Ti^ say, via an expression of the form 

\ rf p = -v m j?„ = n; p + . (so) 

It can easily be seen from (0) that the required deformation tensor will be 
given simply by 

= V p „ \ rf v = - J-W, • (31) 

The middle and last indices of this tensor will evidently have the respective 
properties of tangentiality and orthogonality that are expressible as 

±W CT = 0, H^ a Tf v = Q. (32) 

There is no automatic tangentiality or orthogonality property for the first index 
of the deformation tensor (|3l|), which is thus reducible with respect to the 
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tangential and orthogonally lateral projections (||)([3|) to a sum 



— K)fv ~ L^f (33) 
in which such a property is obtained for each of the parts 

= W» . V = - . (34) 

which satisfy the conditions 

V<?» = X„K?v = = Kfa rf v ■ (35) 

and 

= W = o = ■ ( 36 ) 

The first of these decomposed parts, which might be described as the tan- 
gential turning tensor, is expressible explicitly as 

Kf v = rf a ri\V T rf v , (37) 

and is thus evidently identifiable with the second fundamental tensor as de- 
fined by ([Ti]). The corresponding explicit formula for the second part of the 
decomposition (|3^) is 

L? u = ±^V T J^ . (38) 

Unlike K", which by the Weingartcn identity ( [15"| ) is automatically sym- 
metric on its first two (tangential) indices, the lateral turning tensor, L pl f will in 
general have an antisymmetric part, uj^J 1 say, as well as a symmetric part, 9^ 
say, with respect to its first two (orthogonal) indices. It will be decomposible 
in terms of these parts in the form 

V=V + V' (39) 

with 

V) = ' Vi =0 - ( 4 °) 

Just as the Wiengarten symmetry property (|lj) of K £ is interpretable as 
the integrability condition for the d-surfaces under consideration, so by analogy, 
it is evident that the special condition that is necessary and sufficient for the 
existence of a complementary orthogonal foliation by (n-d)-surfaces is that L^f 
should also be symmetric, i.e. that the rotation tensor uj^jf should vanish. 

The part of L^p that remains even if the foliation is (n-d) surface orthogonal 
is the symmetric part, 9^ , which is the natural generalisation of the usual 
two index divergence tensor 6^ v of an ordinary fluid flow. For a 1-dimcnsional 
timclikc foliation, which will have a unique future directed unit tangent vector 
it^, the contraction 0^ v = 9 pl fu p gives the usual divergence tensor, whose trace 
6 J = V M u M is the ordinary scalar divergence of the flow. The evolution of this 
divergence scalar is governed by the Raychaudhuri identity whose generalisation 
to an evolution equation for 9 pl , was given by Hawking and Ellis ||15| . In the next 
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subsection it will be shown how (following Capovilla and Guven [jl7|) one can 
obtain a further extension J74| ) of the Raychaudhuri identity that provides an 
evolution equation for the three index generalised divergence tensor 9 p which 
(unlike the ordinary divergence tensor 6* M „) is always well defined whatever the 
dimension of the foliation. 

Just as the first order derivatives of the fundamental tensor rf v of a foliation 
are specifiable by the first deformation tensor "H" defined in the previous 
section, so analogously the derivatives of the next order will be fully determined 
by the further specification of a corresponding second deformation tensor, X\ pp 
say, via an expression of the slightly less simple form 

V A H£ P = X ^ P + n \v' H l 7p ~ ^Xp^il'a l ( 41 ) 

that is easily derivable from the definition 

V^WV- ( 42 ) 

using the projection properties (|32|). 

Just as the second fundamental tensor of an individual surface was obtained 
from the first deformation tensor of the foliation by tangential projection accord- 
ing to (|33"|), so analogously the corresponding third fundamental tensor, u , 
as defined by (|24|), will be obtainable from the second deformation tensor by 
the slightly less simple tangential projection operation 

s = v\ < ( V P + h; p ) . (43) 

It can be seen to follow from the Weingarten integrability property ( |I5| ) that 
the antisymmetric part of this third fundamental tensor will be given directly 
by the simple projection operation 

S M v p = V a xV^X [af3 ] P - (44) 

It is easy to use the projection properties (|3^ ) to verify that the tensor obtained 
by this procedure does indeed satisfy the original definition (|i|) of the third 
fundamental tensor, which can thereby be seen to be unambigously determined, 
like the second fundamental tensor, just by the specification of an individual 
imbedded surface, independently of any extension to a space covering foliation. 

The higher order differential information contained in the second deforma- 
tion tensor and the third fundamental tensor is not needed for the evaluation of 
the curvature tensors discussed in the following sections, even though they are 
dependent on the second derivatives of the background metric itself. However 
a useful biproduct of the curvature analysis that follows is the derivation of of 
a separation identity (interpretable as a generalisation of the historic Codazzi 
identity) that expresses the antisymmetric part of the second deformation ten- 
sor - and hence by (|44|) that of the third fundamental tensor - directly as the 
corresponding projection of the background spacetime curvature Tt^u p- 
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1.7 *The adapted connection and curvature of a foliation. 



Due to existence of the decomposition whereby a background spacetime vector, 
with components £ M say, is split up as the sum of its surface tangential part rf- v 
and its surface orthogonal part J^", there will be a corresponding adaptation 
of the ordinary concept of parallel propagation with respect to the background 
connection r„" '. The principle of the adapted propagation concept is to follow 
up an ordinary operation of infinitesimal parallel propagation by the projection 
adjustment that is needed to ensure that purely tangential vectors propagate 
onto purely tangential vectors while purely orthogonal vectors propagate onto 
purely orthogonal vectors. Thus for purely tangential vectors, the effect of the 
adapted propagation is equivalent to that of ordinary internal parallel propaga- 
tion with respect to the induced metric in the imbedded surface, while for purely 
orthogonal vectors it is interpretable as the natural generalisation of the stan- 
dard concept of Fermi- Walker propagation. For an infinitesimal displacement 
dx^ the deviation between the actual component variation (dx v )d y ^ and the 
variation that would be obtained by the corresponding adapted propagation law 
will be expressible in the form [dx v yD v ^ where T> denotes the corresponding 
adapted differentiation operator, whose effect will evidently be given by 

V v e = rf p % ) + ) • (45) 

It can thus be seen that this operation will be expressible in the form 

T>„ C = V„f " + = d v r + A»X , (46) 

where the adapted foliation connection components A.J' ' are given by the for- 
mula 

p = ^p p a n p ' (^U 

in which the ol" are the components of the relevant adaptation tensor, whose 
components can be seen from (^]) to be given by 

a ; p = (^-i;)v^V (48) 

The entirely covariant index lowered version of this adaptation tensor can thus 
be seen to be expressible in terms of the deformation tensor defined by ( |3l| ) in 
the form 

a^p = 2H^[ up ] . (49) 

The fact that this last expression is manifestly antisymmetric with respect to 
the last two indices of the adaptation tensor makes it evident that, like the usual 
Riemannian differentiation operator V, the adapted differentiation operator V 
will commute with index raising or lowering, since the metric remains invariant 
under adapted propagation: 

2^p = 0. (50) 

However, unlike V, the adapted differentiation operator has the convenient prop- 
erty of also commuting with tangential and orthogonal projection, since it can be 
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seen to follow from (^0|) and ( p2| ) that the corresponding operators also remain 
invariant under adapted propagation: 

T>»rT p = 0, V, = 0. (51) 

There is of course a price to be paid in order to obtain this considerable 
advantage of T> over V, but it is not cxhorbitant: all that has to be sacrificed is 
the analogue of the symmetry property 

Vp1=°- (52) 

expressing the absence of torsion in the Riemannian case. For the adapted 
connection A,, l/ n , the torsion tensor defined by 

will not in general be zero. 

When one is dealing not with an extended foliation but just with an individ- 
ual imbedded surface, the general adapted differentiation operator T> will not 
be well defined, but adapted differentiation in tangential directions will still be 
specified by the corresponding unambiguously defined restricted differentiation 
operator T> whose tensorial component representation is given by the adapted 
analogue of (^), namely 

V p = rfpVu . (54) 

Restriction to a single imbedded surface will not invalidate the properties of 
commuting with index raising and tangential or orthogonal projection as ex- 
pressed in the general case by ([SO]) and (|5l|), so in abstract symbolical notation 
we shall have not only [P,g] = 0, but also [2?,??] =0, and \V ,_lJ = 0. It is 



evident from ( [16[ ) and (53) that the effect of this tangcntially restricted adapted 
differentiation operator on any surface supported (but not necessarily surface 
tangential) vector field with components £ M will be given by 

^ r = %e + %%e , < p = 7?>/ p , (55) 

where the relevant tangentially restricted part of the adaptation tensor, which 
may naturally described as the acceleration tensor, is given explicitly as function 
only of the ordinary (restricted) second fundamental tensor of the imbedded 
surface by the simple formula 

a pvp = . (56) 

As a simple illustration, it is to be remarked that in the particular case for 
which the imbedded surface is a particle world-line (with d=l) having a timclikc 
tangent vector ue^ satisfying the unit normalisation condition u v u v = — 1, for 
which the first and second fundamental tensors will be given respectively by 
n^ u = —u^u" and K" = UpU v u p where u p = u a \7 a u p are the components of 
the ordinary acceleration vector, it follows that the corresponding acceleration 
tensor will be given by a pvp = 2u p ,U[ v ii p \. 
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The curvature associated with the adapted connection (|4^) can be read out 
from the ensuing commutator formula, which, for an abitrary vector field with 
components will take the standard form 

v^v v] ^ = - e/„xu' , (57) 

in which the torsion tensor components Qf v are as defined by ([33]) while the 
components T ^ l f a are defined by a Yang-Mills type curvature formula of the 
form 

T^" a = 2d { ,A u] p a + • (58) 

Although the connection components p from which it is constructed are 
not of tensorial type, the resulting curvature components are of course strictly 
tensorial. This is made evident by evaluating the components ( |5S| ) of this amal- 
gamated foliation curvature in terms of the background curvature tensor 

Tl^ a = 2 VV* + T ^J a - , (59) 

and the adaptation tensor a' ' given by ( p9"| ) , which gives the manifestly tenso- 
rial expression 

V. = + 2 + ^ T a v]Ta . (60) 

Although it does not share the full set of symmetries of the Riemann tensor, 
the foliation curvature obtained in this way will evidently be antisymmetric in 
both its first and last pairs of indices: 

J~[lVp(7 ^~\fj,v\\fi&\ ■ (61) 

Using the formula fljsj), it can be see from ( |32| ) and (fhlj), that the difference be- 
tween this adapted curvature and the ordinary background Riemann curvature 
will be given by 

V* - V = 4-V] [H + Mi^K] r + ™ { ;>"H V] Z , (62) 

which superficially appears to depend on the higher order derivatives involved 
in the second deformation tensor X J" 7 . 

This appearance is however deceptive. The reason for qualifying the foliation 
curvature tensor ([58]) as "amalgamated" is that because the adapted derivation 
operator has been constructed in such a way as to map purely tangential vector 
fields onto purely tangential vector fields, and purely orthogonal vector fields 
onto purely orthogonal vector fields, it follows that the same applies to the cor- 
responding curvature, which will therefore consist of an additive amalgamation 
of two separate parts having the form 

•^ilifa = 'Pfj.ifa + Q^ifa > (63) 

in which the "inner" curvature acting on purely tangential vectors is given by a 
doubly tangential projection as 

VV=Va<^, (64) 
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while the "outer" curvature acting on purely orthogonal vectors is given by a 
doubly othogonal projection as 

%u P a = VAA, ■ (65) 

It is implicit in the separation expressed by (^3|) that the mixed tangential and 
orthogonal projection of the adapted curvature must vanish: 

•V\<-4 = 0. (66) 

To get back, from the extended foliation curvature tensors that have just 
been introduced, to their antecedent analogues JlO[ for an individual embedded 
surface, the first step is to construct the amalgamated embedding curvature ten- 
sor, F ' p ' a say, which will be obtainable from the corresponding amalgamated 
foliation curvature T ^ v p a by a doubly tangential projection having the form 

F^ P „=v\M P ^J a - (67) 

As did the extended foliation curvature, so also this amalgated embedding cur- 
vature will separate as the sum of "inner" and "outer" parts in the form 

v. = v. + V* > ( 68 ) 

in which the "inner" embedding curvature is given, consistently wit ( p8| ) by 
another doubly tangential projection as 

R M Sa = VA VTk n\ = V^rja . (69) 



while the "outer" embedding curvature is given, consistently with (|29|) by the 
corresponding doubly lateral projection as 

V. = ^Va4X = v^viQj* ■ (70) 

The non-trival separation identity ( |66| ) can be considered as a generalisation 
to the case of foliations of the relation that is itself interpretable as an extended 
generalisation to higher dimensions of the historic Codazzi equation that was 
originally formulated in the restricted context of 3-dimcnsional flat space. It 
can be seen from (^||) that this extended Codazzi identity is expressible as 

2X btv] P t , + K lo ,\Tf l £, = 0, (71) 

which shows that the relevant higher derivatives are all determined entirely 
by the Riemannian background curvature so that no specific knowledge of the 
second deformation tensor is needed. The formula (^4) can also be used to 
evaluate the "inner" tangential part of the foliation curvature tensor as 

V, = 2?V T ?W + K^ v \ rf K V \ (72) 
and to evaluate the "outer" orthogonal part of the foliation curvature tensor as 

Q,Jc - 2?V P ^]- J^A ■ (73) 
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The formula ([72h for the "inner" foliation curvature is evidently classifiable 



as an extension of the proceeding generalisation ( 28 ) of the historic Gauss equa- 
tion, while the formula ([73) for the "outer" foliation curvature is an extension 
the more recently derived 10 generalisation (|2^) of what has sometimes been 
referred to as the "Ricci equation" but what would seem more appropriately 
describablc as the Schouten equation, with reference to the earliest relevant 
source with which I am familiar G3 . since long after the time of Ricci it was 
not yet understood even by such a leading geometer as Eisenhart [Q . In much 
the same way, the non-trivial separation identity ( |66|) provides the relation ( f7l| ) 
that can be considered as a generalisation to the case of foliations of the identity 
( p5|) that is itself interpretable as an extended generalisation to higher dimen- 
sions of the historic Codazzi equation that was originally formulated in the 
restricted context of 3-dimensional flat space. The corresponding doubly lateral 
projection of the foliation curvature would provide the analogous result for the 
orthogonal foliation by (n-d)-surfaces that would exist in the irrotational case 
for which the lateral turning tensor given by ( |39[ ) is symmetric. Finally the 
corresponding mixed tangential and lateral projection of ( |62| ) gives an identity 
that is expressible in terms of foliation adapted differentiation ([TJ]) as 

+ »W V = K ^\ K ^ + L ^ L x« + ^X/a V p A • (74) 

This last result is interpretable as the translation into the pure background 
tensorial formalism used here of the recently derived generalisation p7[ to higher 
dimensional foliations of the well known Raychaudhuri equation (whose original 
scalar version, and its tensorial extension [ pL5| , were formulated just for the 
special case of a foliation by 1-dimensional curves). The complete identity (^TJ) 
is therefore interpretable as an amalgamated Raychaudhuri-Codazzi identity. 



1.8 *The special case of a string worldsheet in 4-dimensions 

The application with which we shall mainly be concerned in the following work 
will be the case d=2 of a string. An orthonormal tangent frame will consist 
in this case just of a timclikc unit vector, lq^, and a spacelike unit vector, 
ti^j whose exterior product vector is the frame independent antisymmetric unit 
surface element tensor 

£»» = 2i^l^ = 2(H 7 ir 1/2 ^ < ] 1 > (75) 

whose tangential gradient satisfies 

V\£ ^ = -2K Xp ^E v]p . (76) 

(This is the special d=2 case of what is, as far as I am aware, the only wrongly 
printed formula in the more complete analysis [flo| on which this section is based: 
the relevant general formula (B9) is valid as printed only for odd d, but needs 
insertion of a missing sign adjustment factor (— l) d_1 in order to hold for all 
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d.) In this case the inner rotation pseudo tensor jlO| ) is determined just by a 
corresponding rotation covector p M according to the specification 

PS* = k £ "»P*> Px = P^£\- (77) 

This can be used to see from ( [j"2| ) that the Ricci scalar, 

R = R? v R^v — R pll f '„ , (78) 

of the 2-dimensional worldshcet will have the well known property of being a 
pure surface divergence, albeit of a frame gauge dependent quantity: 

R = %(E» v p u ). (79) 

In the specially important case of a string in ordinary 4-dimensional spacctimc, 
i.e. when we have not only d=2 but also n=4, the antisymmetric background 
measure tensor £ A ^p can b e use d to determine a scalar (or more strictly, since 
its sign is orientation dependent, a pseudo scalar) magnitude fl for the outer 
curvature tensor (jl^) (despite the fact that its traces are identically zero) ac- 
cording to the specification 

Q = \n Xilvp E^ v P . (80) 

Under these circumstances one can also define a "twist" covector zu^, that is 
the outer analogue of p M , according to the specification 

w v = | w/ A e Awot £P° . (81) 



This can be used to deduce from (|13|) that the outer curvature (pseudo) scalar 
Q of a string worldshcet in 4-dimcnsions has a divergence property of the same 



kind as that of its more widely known Ricci analogue (79), the corresponding 
formula being given by 

(2=%(£" v m v ). (82) 



It is to be remarked that for a compact spacclike 2-surface the integral of ( |76| ) 
gives the well known Gauss Bonnet invariant, but that the timelike string world- 
sheets under consideration here will not be characterised by any such global 
invariant since they will not be compact (being open in the time direction even 
for a loop that is closed in the spacial sense). The outer analogue of the Gauss 
Bonnet invariant that arises from (^o|) for a spacelike 2-surface has been dis- 
cussed by Penrose and Rindlcr Jl£| but again there is no corresponding global 
invariant in the necessarily non-compact timelike case of a string worldshcet. 

2 Laws of motion for a regular pure brane com- 
plex 
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2.1 Regular and distributional formulations of a brane ac- 
tion 

The term p-brane has come into use B, jl9| to describe a dynamic system 
localised on a timelike support surface of dimension d=p+l, imbedded in a 
spacetime background of dimension n>p. Thus at the low dimensional extreme 
one has the example of a zero - brane, meaning what is commonly referred 
to as a "point particle" , and of a 1-brane meaning what is commonly referred 
to as a "string" . At the high dimensional extreme one has the "improper" 
case of an (n-l)-brane, meaning what is commonly referred to as a "medium" 
(as exemplified by a simple fluid), and of an (n-2)-brane, meaning what is 
commonly referred to as a "membrane" (from which the generic term "brane" 
is derived). A membrane (as exemplified by a cosmological domain wall) has 
the special feature of being supported by a hypcrsurface, and so being able to 
form a boundary between separate background space time regions; this means 
that a 2-brane has the status of being a membrane in ordinary 4-dimcnsional 
spacetime (with n = 4) but not in a higher dimensional (e.g. Kaluza Klein type) 
background. 

The purpose of the present section is to consider the dynamics not just of 
an individual brane but of a brane complex or "rigging model" (^] such as is 
illustrated by the nautical archetype in which the wind - a 3-brane - acts on 
a boat's sail - a 2-branc - that is held in place by cords - 1-branes - which 
meet at knots, shackles and pulley blocks that are macroscopically describable 
as point particles - i.e. 0-branes. In order for a a set of branes of diverse 
dimensions to qualify as a "geometrically regular" brane complex or "rigging 
system" it is required not only that the support surface of each (d-l)-brane 
should be a smoothly imbedded d-dimensional timelike hypersurface but also 
that its boundary, if any, should consist of a disjoint union of support surfaces of 
an attatched subset of lower dimensional branes of the complex. (For example 
in order qualify as part of a regular brane complex the edge of a boat's sail can 
not be allowed to flap freely but must be attatched to a hem cord belonging 
to the complex.) For the brane complex to qualify as regular in the strong 
dynamic sense that will be postulated in the present work, it is also required 
that a member p-brane can exert a direct force only on an an attached (p— 1)- 
brane on its boundary or on an attached (p+l)-brane on whose boundary it is 
itself located, though it may be passively subject to forces exerted by a higher 
dimensional background field. For instance the Peccei-Quin axion model gives 
rise to field configurations rcprcscntable as regular complexes of domain walls 
attached to strings JtJ, 20 , and a bounded (topological or other) Higgs 



vortex defect terminated by a pair of pole defects |22j , |23| , J24J , |2jJ , |26| , 
may be represented as a regular brane complex consisting of a finite cosmic 
string with a pair of point particles at its ends, in an approximation neglecting 
Higgs field radiation. (However allowance for radiation would require the use 
of an extended complex including the Higgs medium whose interaction with 
the string - and a fortiori with the terminating particles - would violate the 
regularity condition: the ensuing singularities in the back reaction would need 
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to be treated by a renormalisation procedure of a kind pl| , p8[ , p9[ , j30| 
whose development so far has been beset with difficulties in preserving exact 
local Lorentz invariance, an awkward problem that is beyond the scope of the 
present article.) 

The present section will be restricted to the case of a brane complex that is 
not only regular in the sense of the preceeding paragraph but that is also pure 
(or "fine" ) in the sense that the lengthscalcs characterising the internal structure 
of the (defect or other) localised phenomenon represented by the brane models 
are short compared with those characterising the macroscopic variations under 
consideration so that polarisation effects play no role. For instance in the case 
of a point particle, the restriction that it should be describable as a "pure" zero 
brane simply means that it can be represented as a simple monopole without 
any dipole or higher multipole effects. In the case of a cosmic string the use 
of a "pure" 1-brane description requires that the underlying vortex defect be 
sufficiently thin compared not only compared with its total length but also 
compared with the lengthscalcs characterising its curvature and the gradients 
of any currents it may be carrying. The effect of the simplest kind of curvature 



corrections beyond this "pure brane" limit will be discussed in Subsection 3.5 
but in the rest of these lectures, as in the present section, it will be assumed 
that the ratio of microscopic to macroscopic lengthscales is sufficiently small for 
description in terms of "pure" p-branes to be adequate. 

The present section will not be concerned with the specific details of par- 
ticular cases but with the generally valid laws that can be derived as Noether 
identities from the postulate that the model is governed by dynamical laws 
derivable from a variational principle specified in terms of an action function 
I. It is however to be emphasised that the validity at a macroscopic level of 
the laws given here is not restricted to cases represented by macroscopic models 
of the strictly conservative type directly governed by a macroscopic variational 
principle. The laws obtained here will also be applicable to classical models of 
dissipative type (e.g. allowing for resistivity to relative flow by internal currents) 
as necessary conditions for the existence of an underlying variational description 
of the microscopic (quantum) degrees of freedom that are allowed for merely as 
entropy in the macroscopically averaged classical description. 

In the case of a brane complex, the total action T will be given as a sum 
of distinct d-surface integrals respectively contributed by the various (d-1)- 
branes of the complex, of which each is supposed to have its own corresponding 
Lagrangian surface density scalar < d '£ say. Each supporting d-surface will be 
specified by a mapping a i— > x{a} giving the local background coordinates x^ 
(/i=0, .... , n-1) as functions of local internal coordinates a % ( i=0, ... , d-1). 
The corresponding d-dimcnsional surface metric tensor ( d >7y that is induced (in 



the manner described in Subsection 1.1) as the pull back of the n-dimensional 



background spacetime metric g^ vi will determine the natural surface measure, 
< d >dS, in terms of which the total action will be expressible in the form 



E 



J < d >dS < d >£ , < d >dS = VTI^tII d d a . (83) 
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As a formal artifice whose use is an unnecessary complication in ordinary dy- 
namical calculations but that can be useful for purposes such as the calculation 
of radiation, the confined (d-surface supported) but locally regular Lagrangian 
scalar fields < d >£ can be replaced by corresponding unconfined, so no longer regu- 
lar but distributional fields < d >£, in order to allow the the basic multidimensional 
action (^3|) to be represented as a single integral, 

J=/dS^< d >£, dS = y^\\gj\d n x. (84) 

^ d 

over the n-dimcnsional background spacctime. In order to do this, it is evi- 
dent that for each (d-l)-brane of the complex the required distributional action 
contribution < d >£ must be constructed in terms of the corresponding regular 
d-surface density scalar < d >£ according to the prescription that is expressible in 
standard Dirac notation as 

< d >£ = ||g|r 1/2 J < d >dS ^CS n [x - x{a}} . (85) 

2.2 Current, vorticity, and stress-energy tensor 

In the kind of model under consideration, each supporting d-surface is supposed 
to be endowed with its own independent internal field variables which are al- 
lowed to couple with each other and with their derivatives in the corresponding 
d-surfacc Lagrangian contribution < d >£, and which are also allowed to couple 
into the Lagrangian contribution ( d - 1 >£ on any of its attached boundary (d-1) 
surfaces, though - in order not to violate the strong dynamic regularity con- 
dition - they are not allowed to couple into contributions of dimension (d-2) 
or lower. As well as involving its own d-brane surface fields and those of any 
(d+f ) brane to whose boundary it may belong, each contribution < d >£ may also 
depend passively on the fields of a fixed higher dimensional background. Such 
fields will of course always include the background spacctime metric g pv itself. 
Apart from that, the most commonly relevant kind of backround field (the only 
one allowed for in the earlier analysis, ||) is a Maxwellian gauge potential 
whose exterior derivative is the automatically "closed" electromagnetic field, 

F»v = 2^.4,] , V[»F up] = . (86) 

Although many other possibilities can in principle be enviseaged, for the sake 
of simplicity the following analysis will not go beyond allowance for the only 
one that is important in applications to the kind of cosmic or supcrfluid defects 
that are the subject of discussion in the present volume, namely an antisym- 
metric Kalb-Ramond gauge field B^ v = —B VI1 whose exterior derivative is an 
automativally closed physical current 3-form, 

N^ P = ^B vp] , V^N upa] = . (87) 

Just as a Maxwellian gauge transformation of the form A p \— ► A p + V^a for an 
arbitrary scalar a leaves the electromagnetic field (p6|) invariant, so analogously 
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a Kalb-Ramond gauge transformation B ^ i— ► B ^ + 2V[ M Xi/] fo r an arbitrary 
covector \n leaves the corresponding current 3-form ( |S7| ) invariant. In applica- 
tions to ordinary 4-dimensional spacetime, the current 3-form will just be the 
dual N^vp — EpupaN of an ordinary current vector iV M satisfying a conserva- 
tion law of the usual type, V M iV M = 0. Such a Kalb-Ramond representation can 
be used to provide an elegant variational formulation for ordinary perfect fluid 
theory plj and is particularly convenient for setting up "global" string models 
of vortices both in a simple cosmic axion or Higgs field fffi] , J33|, J34| and in a 
superfluid |55| such as liquid Hclium-4. 

In accordance with the preceeding considerations, the analysis that follows 
will be based on the postulate that the action is covariantly and gauge invari- 
antly determined by specifying each scalar Lagrangian contribution < d )£ as a 
function just of the background fields, A M , B^ v and of course g^, and of any 
relevant internal fields (which in the simplest non-trivial case - exemplified by 
string models J|6| , J|lj of the category needed for the macroscopic description 
of Wittcn type [0 superconducting vortices - consist just of a phase scalar ip). 
In accordance with the restriction that the branes be "pure" or "fine" in the 
sense explained above, it is postulated that polarisation effects are excluded by 
ruling out couplings involving gradients of the background fields. This means 
that the effect of making arbitrary infinitesimal "Lagrangian" variations 5Ap, 
SB^v , 8gp U of the background fields will be to induce a corresponding variation 
IT of the action that simply has the form 

SI = Y, [ (d) ^ { w ~f5-^ + \ (dy w^5Bp„ + | < d 'T^^^} , 
d J 

provided either that that the relevant independent internal field components are 
fixed or else that the internal dynamic equations of motion are satisfied in accor- 
dance with the variational principle stipulating that variations of the relevant 
independent field variables should make no difference. For each d-brane of the 
complex, this partial differentiation formula implicitly specifies the correspond- 
ing electromagnetic surface current density vector <d) j M , the surface vorticity 
flux bivector Ww^" — —^w^^, and the surface stress momentum energy den- 
sity tensor WT^" = WT ufl , which are formally expressible more explicitly as 

^=ax;> (d ^ = 2 ^w (89) 



and 



7 ' = 2- + W£C<V U , (90) 

og^ 



of which the latter is obtained using the formula 

5{^dS)= 1 -^{8gp_ v )^dS, (91) 

where i^qf 1 " is the rank-d fundamental tensor of the d-dimensional imbedding, 
as defined in the manner described in the preceeding section. 
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2.3 Conservation of current and vorticity 

The condition that the action be gauge invariant means that if one simply sets 
SA^ = V^a, SB^y = 2V[ M x^]i d 9^ = 0, for arbitrarily chosen a and x^ then 
61 should simply vanish, i.e. 

E fd^S { (d >7 M V M a + t*Wi""VnXv} = • (92) 

d J 

In order for this to be able to hold for all possible fields a and Xn it is evident that 
the surface current < d) j M and the vorticity flux bivcctor must (as one would 

anyway expect from the consideration that they depend just on the relevant 
internal d-surface fields) be purely d-surface tangential, i.e. their contractions 
with the relevant rank (n— d) orthogonal projector < d l/^ = g^ v — Wr/^ must 
vanish: 

(d )L^ W f = o 5 tdjLM <d)^p = o . (93) 

Hence, decomposing the full gradient operator V M as the sum of its tangentially 
projected part < d 'V M = ^rf^v and of its orthogonally projected part tgl^y 
and noting that by (|93| ) the latter will give no contribution, one sees that ( p2| ) 
will take the form 

E / (d) ^ { (d % ( (d) F" + v^xi) - a {d) %f - xu (d % ^r"} = , 

d ^ 

| (94) 

in which first term of each integrand is a pure surface divergence. Such a 
divergence can be dealt with using Green's theorem, according to which, for 
any d-dimcnsional support surface (d \S of a (d-l)-brane, one has the identity 

J W(6( d % {d Y = <j> ^-^dS^Xf, , (95) 

where the integral on the right is taken over the boundary (d-l)-surface of d < d \S 
of ( d \S, and < d 'A M is the (uniquely defined) unit tangent vector on the d-surface 
that is directed normally outwards at its (d-l)-dimensional boundary. Bearing 
in mind that a membrane support hypersurface can belong to the boundary of 
two distinct media, and that for d< n-3 a d-branc may belong to a common 
boundary joining three or more distinct (d+l)-brancs of the complex under 
consideration, one sees that (B4) is equivalent to the condition 



^2 (v)dS (a(wV li Mj' t - E ( % (d) i' 1> 

P "* d=p+l 

+X»( (p) % (py w fny - E (d ' A ^ (dmjfU/ ) } = °' ( 96 ) 



d=p+l 

where, for a particular p-dimcnsionally supported (p-l)-brane, the summation 
"over d=p+l" is to be understood as consisting of a contribution from each 
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(p+l)-dimcnsionally supported p-brane attached to it, where for each such p- 
brane, ( d >A M denotes the (uniquely defined) unit tangent vector on its (p+1)- 
dimcnsional support surface that is directed normally towards the p-dimensional 
support surface of the boundary (p-l)-branc. The Maxwell gauge invariance 
requirement to the effect that (^6|) should hold for arbitrary a can be seen to 
entail an electromagnetic charge conservation law of the form 

<p % = (d) -V <d) F ■ ( 97 ) 

d=p+l 

This can be seen from ( ]95| ) to be be interpretable as meaning that the total 
charge flowing out of particular (d-l)-brane from its boundary is balanced by 
the total charge flowing into it from any d-branes to which it may be attached. 
The analogous Kalb-Ramond gauge invariance requirement that ( |96| ) should also 
hold for arbitrary \n can be seen to entail a corresponding vorticity conservation 
law of the form 

w% = Y (d)/ V (d)w ' 11/ ■ ( 98 ) 

d=p+l 

A more sophisticated but less practical way of deriving the foregoing conserva- 
tion laws would be to work not from the expression ( |S3| ) in terms of ordinary 
surface integrals but instead to use the superficially simpler expression (Si]) in 



terms of distributions, which leads to the replacement of (97) by the ultimately 



equivalent (more formally obvious but less directly meaningful) expression 

'Wji*) = (99) 



involving the no longer regular but Dirac distributional current WjP that is 
given in terms of the corresponding regular surface current by 

<d) .f' = ll.9ir 1/2 / (d) ^ (d) F<S n [x - x{a}} . (100) 



Similarly one can if one wishes rewrite the vorticity flux conservation law ( pq ) 
in the distributional form 



=0, (101) 



where the distributional vorticity flux Ww 11 " is given in terms of the correspond- 
ing regular surface flux ( d W fJ " L ' by 

< d >u^ = ||.g|r 1/2 / < d >dS ^w^S n [x - x{a}} . (102) 



It is left as an entirely optional exercise for any readers who may be adept in 
distribution theory to show how the ordinary functional relationships (97) and 
(38h can be recovered by by integrating out the Dirac distributions in (99) and 



(101) 
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2.4 Force and the stress balance equation 

The condition that the hypothetical variations introduced in ( |S8| ) should be 
"Lagrangian" simply means that they are to be understood to be measured 
with respect to a reference system that is comoving with the various brancs 
under consideration, so that their localisation with respect to it remains fixed. 
This condition is necessary for the variation to be meaningly definable at all for 
a field whose support is confined to a particular brane locus, but in the case of 
an unrestricted background field one can enviseage the alternative possibility 
of an "Eulerian" variation, meaning one defined with respect to a reference 
system that is fixed in advance, independently of the localisation of the brane 
complex, the standard example being that of a Minkowski reference system 
in the case of a background that is flat. In such a case the relation between 
the more generally meaningfull Lagrangian (comoving) variation, denoted by 5, 
and the corresponding Eulerian (fixed point) variation denoted by S say will be 
given by Lie differentiation with respect to the vector field £ M say that specifics 
the infinitesimal of the comoving reference system with respect to the fixed 
background, i.e. one has 

S-6= ££, (103) 

where the Lie differentiation operator ££ is given for the background fields 
under consideration here by 

££A„ = rv,.4 M + A v Vnt v , (104) 

{£B^ = + 2B p[y V^ p , (105) 

= 2V (Al ^) . (106) 

This brings us to the main point of this section which is the derivation of the 
dynamic equations governing the extrinsic motion of the branes of the complex, 
which are obtained from the variational principle to the effect that the action X 
is left invariant not only by infinitesimal variations of the relevant independent 
intrinsic fields on the support surfaces but also by infinitesimal displacements 
of the support surfaces themselves. Since the background fields A^, £? M „, and 
g^u are to be considered as fixed, the relevant Eulerian variations simply vanish, 
and so the resulting Lagrangian variations will be directly identifiable with the 



corresponding Lie derivatives - as given by (106) - with respect to the gener- 
ating vector field £ M of the infinitesimal displacement under consideration. The 
variational principle governing the equations of extrinsic motion is thus obtained 
by setting to zero the result of substituting these Lie derivatives in place of the 
corresponding Lagrangian variations in the more general variation formula ( |S8| ) , 
which gives 

f (d) <£ { W T + i ^ v {£B^ + \ < d >T^ feg^} = . (107) 

The requirement that this should hold for any choice of £ M evidently implies that 
the tangentiality conditions for the surface fluxes ^j* 1 and ( d >?U A " y must be 
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supplemented by an analogous d-surface tangentiality condition for the surface 
stress momentum energy tensor ( d )T MI/ , which must satisfy 

CdlM WT" P = . (108) 

(as again one would expect anyway from the consider ation that it depends just 
on the relevant internal d-surface fields). This allows (106) to be written out in 
the form 

Y, I w &{e(F pll \N ptlv 

- < d % ' d 'T^ - A p < d % < d )J M - B pv < d % < d >uy^) 

+ (d %(e(A p ^j>*+B p „Ww^+ < d >T%))} = 0, (109) 

in which the final contribution is a pure surface divergence that can be dealt with 
using Green's theorem as before. Using the results ( |97|) and ( ^8|) of the analysis 
of the consequences of gauge invarianc e an d proceeding as in their derivation 
above, one sees that the condition for (|l09| ) to hold for an arbitrary field £ M is 
that, on each (p-l)-brane of the complex, the dynamical equations 

M%b>W t p =Wf p , (110) 

should be satisfied for a total force density (p '/ p given by 

(p) f P = (p) 7p + M f P , (in) 

where iv) f p is the contribution of the contact force exerted on the p-surface by 
other members of the brane complex, which takes the form 

(p) /p= X! w X M < d *T» P , (112) 

d=p+l 

while the other force density contribution Cp)/ p represents the effect of the ex- 
ternal background fields, which is given by 

(p, 7p = F pp <p>J" + ±N pia , ^w^ . (113) 



As before, the summation "over d=p+l" in (112) is to be understood as con- 
sisting of a contribution from each of the p-branes attached to the (p-1)- brane 
under consideration, where for each such attached p-brane, ( d 'Au denotes the 
(uniquely defined) unit tangent vector on its (p+l)-dimensional support sur- 
face that is directed normally towards the p-dimensional support surface of the 
boundary (p-l)-brane. The first of the background force contibutions in (113) is 
of course the Lorentz type force density resulting from the effect of the electro- 
magnetic field on the surface current, while the other contribution in ( p_13[ ) is a 
Joukowsky type force density (of the kind responsible for the lift on an aerofoil) 
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resulting from the Magnus effect, which acts in the case of a "global" string p2j, 
[ p3| though not in the case of a string of the "local" type for which the relevant 
vorticity flux <p)zU M!/ will be zero. As with the conservation laws (|97]) and (|9^), so 
also the explicit force density balance law expressed by ( |l 10| ) can alternatively 
be expressed in terms of the corresponding Dirac distributional stress momen- 
tum energy and background force density tensors, (A)T^ U and ^'/p, which are 
given for each (d-l)-brane in terms of the corresponding regular surface stress 
momentum energy and background force density tensors wT^ u and 



<d) /„ by 



< d >T^ = || 5 !r 1/2 / <d, ^S ^T^5 n [x - x{a}} 



and 



(d) A 



ur 1 ' 2 ^ds^f^[x~x{<j}\ 



(114) 



(115) 



The equivalent - more formally obvious but less explicitly meaningful - distri- 
butional versional version of the force balance law flllOp takes the form 



(116) 



where the total Dirac distributional force density is given in terms of the electro- 
magnetic current distributions ( |100| ) and the vorticity flux distributions (102) 
by 

d d 

It is again left as an optional exercise for readers who are adept in the use of 
Dirac distributions to show that the system ( 1 1C ) , (112) , (|113| ) is obtainable 
from ( |lT6| ) and ( pTfj ) by substituting ( |i~00| ), (^02^7 (P^4^ (P^K 

As an immediate corollary of (11C), it is to be noted that for any vector field 
I' 1 that generates a continuous symmetry of the background spacetime metric, 
i.e. for any solution of the Killing equations 



0, (118) 
one can construct a corresponding surface momentum or energy density current 



that will satisfy 



(p>V M ( 



(Pipe — (PiT^i 



p,pM = (d); v (d,pAI + <p) /^ -Al 

d=p+l 



(119) 



(120) 



In typical applications for which the n-dimensional background spacetime can 
be taken to be flat there will be n independent translation Killing vectors which 
alone (without recourse to the further n(n-l)/2 rotation and boost Killing vec- 
tors of the Lorentz algebra) will provide a set of relations of the form (120) that 
together provide the same information as that in the full force balance equation 



(110) or (116) 
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2.5 The equation of extrinsic motion 



Rather than the distributional version (116), it is the explicit version (110) of 



the force balance law that is directly useful for calculating the dynamic evolution 



of the brane support surfaces. Since the relation (116) involves n independent 
components whereas the support surface involved is only p-dimensional, there 
is a certain redundancy, which results from the fact that if the virtual displace- 
ment field £ M is tangential to the surface in question it cannot affect the action. 



Thus if ( p i^ t/ ^ I/ = 0, the condition (107) will be satisfied as a mere identity 



provided of course that the field equations governing the internal fields of the 
system are satisfied. It follows that the non-rcdundent information governing 
the extrinsic motion of the p-dimensional support surface will be given just by 
the orthogonally projected part of ( |110| ). Integrating by parts, using the fact 
that, by (||) and (17), the surface gradient of the rank-(n-p) orthogonal projec- 



tor Cp)^ will be given in terms of the second fundamental tensor of the 

p-surface by 

(p)V(p1! / =-(p)K p - (p)R p , (121) 
it can be seen that the e xtrinsic equations of motion obtained as the orthogonally 



projected part of (110) will finally be expressible by 

wf^ WKpS = Wt? M <p'/ M . (122) 

It is to be emphasised that the formal validity of the formula that has just 
been derived is not confined to the variational models on which the above 
derivation is based, but also extends to dissipativc models (involving effects 
such as external drag by the background medium |4C|] , |3§| ], |39| ] or mutual re- 
sistance between independent internal currents). The condition that even a 
non-conservative macroscopic model should be compatible with an underlying 
microscopic model of conservative type requires the existence (representing to 
averages of corresponding microscopic quantities) of appropriate stress momen- 



tum energy density and force density fields satisfying (122). 

The ubiquitously applicable formula ( |122| )is intcrprctablc as being just the 
natural higher generalisation of "Newton's law" (equating the product of mass 
with acceleration to the applied force) in the case of a particle. The surface 
stress momentum energy tensor, MT^" , generalises the mass, and the second 
fundamental tensor, ^K^, generalises the acceleration. 

The way this works out in the 1-dimcnsional case of a "pure" point particle 
(i.e. a monopole) of mass to, for which the Lagrangian is given simply by 
<!)£ = —to, is as follows. The 1-dimensional energy tensor will be obtained 
in terms of the unit tangent vector (u^u^ = — 1) as u>7^ ll/ = mu^u^, and 
in this zero-brane case, the first fundamental tensor will simply be given by 
Ci^/jk _ — u ^ u v ^ go that the second fundamental tensor will be ob taine d in 
terms of the acceleration u 11 = u y \/ u u p as Wlf ? = u^u v ii p . Thus (122) can 



be seen to reduce in the case of a particle simply to the usual familiar form 
mil? =W J?. («f. 

(J, J 
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3 Perturbations and curvature effects beyond 
the pure brane limit 

3.1 First order perturbations and the extrinsic character- 
istic equation 

Two of the most useful formulae for the analysis of small perturbations of a string 
or higher brane worldsheet are the expressions for the infinitesimal Lagrangian 
(comoving) variation of the first and second fundamental tensors in terms of 
the corresponding comoving variation Sg pl , of the metric (with respect to the 
comoving reference system). For the first fundamental tensor one easily obtains 

Sf = -V^Sgp* , Sf v - rr±lSg pa (123) 

and, by substituting this in the defining relation (|l4|), the corresponding La- 
grangian variation of the second fundamental tensor is obtained p4| as 

6K p J> = f,\ + {2^K v) T " - K^rf") Sg ar , (124) 

where the Lagrangian variation of the connection ( |20| ) is given by the well known 
formula 

5T a \ = g Xp {V (a 8g T)p - \v p 5g aT ) . (125) 

Since we are concerned here only with cases for which the back ground is fixed in 
advance so that the Eulcrian variation d E will vanish in ( 103 ), the Lagrangian 
variation of the metric will be given just by its Lie derivative with respect to 
the infinitesimal displacement vector field £ M that generates the displacement of 
the worldsheet under consideration, i.e. we shall simply have 

6g aT = 2V ((T £ r) . (126) 



It then follows from (125) that the Lagrangian variation of the connection will 
be given by 

ST a \ = V (ff V T) e A - K\ aT)p ¥ , (127) 

where lZ\ Tp is the background Riemann curvature (which will be negligible 
in typical applications for which the lengthscalcs characterising the geometric 
features of interest will be small compared with those characterising any back- 
ground spacetime curvature) . The Lagrangian variation of the first fundamental 
tensor is thus finally obtained in the form 

5if v = -27 ? ^ l V ,y) r , (128) 
while that of the second fundamental tensor is found to be given by 

SK p f = ^(V^f* - ^f v) n\ Tp i p - A- (Miy) V CT £ A )+ (129) 

{2^RP - gP T K p j) (V,r + V T C.) . (130) 
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It is instructive to apply the forgoing formulae to the case of a free pure 
brane worldsheet, meaning one for which there is no external force contribution 
so that the equation of extinsic motion reduces to the form 



rpfiv TS p 



0. 



(131) 



On varyi ng th e relation (131) using (130) in c onju nction with the orthogonality 
property (10S) and the unperturbed equation ( 131 ) itself, the equation governing 
the propagation of the infinitesimal displacement vector is obtained in the form 



-K PST"" 



(132) 



The extrinsic perturbation equation (132) is by itself only part of the com- 
plete system of perturbation equations governing the evolution of the brane, 
the remaining equations of the system being those governing the evolution of 
whatever surface current pl| and other relevant internal fields on the support- 
ing worldsheet may be relevant. The perturbations of such fields are involved in 



the source term on the right of (132), whose explicit evaluation depends on the 
specific form of the relevant currents or other internal fields. However it is not 
necessary to know the specific form of such internal fields for the purpose just 
of deriving the characteristic velocities of propagation of the extrinsic propaga- 
tions represented by the displacement vector so long as they con tribu te to 
the source term on the right of the linearised perturbation equation (132) only 
at first differential order, so that the characte ristic velocities will be completely 
determined by the first term on the left of ( 132 ) which wil l be the only sec- 
ond differential order contribution. It is apparent from (132) that under these 
conditions the equation for the characteristic tangent covector \n sa Y wn l be 
given independently of any details of the surface currents or other internal fields 
simply §] by 

T^ X „Xu = 0- (133) 

(It can be seen that the unperturbed surface stress momentum energy density 
tensor T^ u plays the same role here as that of the unperturbed metric tensor 
in the analogous characteristic equation for the familiar case of a massless 
background spacetime field, as exemplified by electromagnetic or gravitational 
radiation.) 



3.2 *Higher order displacements by "straight" transporta- 
tion 

Mere linear perturbation theory, such as described in the preceeding subsection, 
is of course insufficient for many physical purposes (such as the treatment of 
gravitational radiation reaction for which it is well known to be necessary to go 
even beyond second order in the post Newtonian approximation) . Allowance for 
actual physical effects involving perturbations at quadratic or higher order is 
beyond the scope of the present notes. Nevertheless even when one's attention 
is restricted to first order in so far as "physical" (on shell) perturbations are 
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concerned, it is still necessary to go to quadratic order in the treatment of 
"virtual" (off shell) perturbations if one wants a variational formulation of the 
linearised dynamical equations. 

The treatment of higher order perturbations - whether "physical" or "vir- 
tual" - is intrinsically straightforcward (though it may of course involve very 
heavy algebra) when it only involves continuous background fields for which 
a purely "Eulerian" (effectively fixed) reference system can be used. However 
when point particles, strings, or higher branes are involved one needs to use 
a comoving "Lagrangian" type reference system, in the manner developed in 
Section ^. In so far as only first order perturbation theory is concerned, the nec- 
essary adjustment between Eulerian and Lagrangian systems can be treated in 



accordance with (103) by the well known procedure whereby the corresponding 
differences are obtained as Lie derivatives with respect to the relevant displace- 
ment vector field. However a new issue arises when one comes to generalise 
such a procedure to quadratic and higher order perturbation theory because 
in a general mathematical context it is only to first order that a small coordi- 
nate displacement will unambiguously characterise and be characterised by a 
corresponding vector. 

Fortunately for what follows, the physical context with which the present 
notes are concerned is not of the most general mathematical kind, but one in 
which the background under consideration is endowed with a preferred connec- 
tion, TJ 1 say, namely the one derived according to the standard Riemannian 
prescription ( p0| ) from the spacetime metric g M „. So long as one has a pre- 
ferred connection (whether Riemannian or not) to work with, there will after 
all be a corresponding precription, namely the one naturally specified by the 
corresponding affinely parametrised geodesic, which we shall refer to simply as 
"straight" transportation, whereby even to second and higher order, a given 
vector will unambiguously determine a corresponding displacement and also (at 
least provided the displacement is not unduly large) vice versa. 

To be explicit, let £^ be the components of a vector at a given position 
with coordinates and let e be the corresponding affinc parameter along 
the geodesic starting from x^ for which the initial affinely normalised tangent 
vector is £ M . This will characterise a well defined displacement x M i— > iL 
specified by setting e = 1 in the solution x j 1 1 with initial conditions 



CM i"\t 4- Vi /"i ri +-H n s~\ \ t t i-» n i-i-i -m /-vJ- *rvt»i ciArl (TA/ 

{0} 



irli = of the affinely parameterised geodesic equation 



£ {e}+ r {eW i {e}*?e} = 0> ( 134 ) 

where a dot denotes the derivative with respect to the affine parameter. The 
required solution of this geodesic equation is conveniently expressible by a Taylor 
expansion whose expression to third order will have the form 

2 3 

xfe - a" = e $x" + e -% 2 x» + ^$ 3 ^ + 0{e 4 } , (135) 

using a systematic notation scheme whereby the symbol $ is introduced for 
the purpose of denoting differentiation with respect to the affine parameter at 
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its initial value, e = 0, so that for example we shall have the abbreviation 
$2;^ = i{ y Proceeding by successive approximations, it is a straightforward 
exercise to verify that, up to third order, the required coefficients will be given 

by 

$x" = e , 8V = -Tf v £ x t x , (136) 

$ 3 *^ = {r^ v) r p \ r^jrrr . (137) 

The purpose of this scheme is to obtain a systematic calculus that can be 
used to obtain the generalisation to higher orders of the standard formula ( [10 3D 
whereby starting from a first order Eulerian variation operator S the corre- 
sponding Lagrangian variation operator <5 as specified with respect to a dis- 
placed reference system will be given by o — 6 = £ where £ is the operator of 
Lie differentiation with respect to the vector £^ that characterises the relative 
displacement. Starting from a finite Eulerian variation operator A allowing for 
higher order corrections obtained by successive applications of the correspond- 
ing infinitesimal differential operator e? L via an Taylor expansion of the form 
A = S + 1 /oS 2 + ... we want to be able to construct the corresponding finite 

EE E 

Lagrangian variation operator A = <5 + 1 /2<5 2 + ... as specified with respect to 
a reference system that has been subject to a relative transportation operation 
representable in terms of local coordinates as the mapping x M i— > ^jiy By ex- 
pressing this mapping in terms of the displacement vector £ M in the manner 
that has just beeen described, the required difference will be obtainable as a 
Taylor expansion in powers of the vector £ M , and of its first (but not higher) 
order gradient tensor V,,^, in the form 

A L -A E =$ + i$ 2 + |$ 3 + ... (138) 

where $, $ 2 , and so on are the relevant operators of first, second, and higher 
order "straight" differentiation with respect to the £ M . The way to work out the 
explicit form of the action of such "straight" differentiation operators on the 
fields of the kind that is of interest - of which the most important is of course 
the metric - will be described in the following subsection. 

Before proceeding, it is to be observed that the geodesic or "straight" trans- 
portation operation, based in the manner just described on the solution of the 



second order equation (134), needs to be distinguished from the corresponding 
Lie transportation operation, which is instead based on the solution of a first or- 
der equation of the form = . It would in principle be possible to develop 
higher order displacement perturbation theory in terms of Lie transportation, 
but the disadvantage of such an aproach is is that (as the price of avoiding de- 
pendence on the preferred connection) it would introduce an inordinately high 
degree of gauge dependence, since it would depend on the specification of the 
vector not just at the initial point but all the way along the trajectory of 
the displacement. However the distinction arises only at quadratic and higher 
order: at first order there is no difference between the "straight" differentiation 
operator denoted by $, that has just been defined and the corresponding Lie 
differentiation operation as customarily denoted by the symbol £. 
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It is to be remarked by the way that, just as is customary to use the more 
explicit notation ££ or alternatively £^ whenever it is necessary to indicate the 
particular vector field £ M with respect to which the Lie differentiation operation 
£ is defined, so analogously we can use the more explicit notation £$ or alter- 
natively $j to indicate the particular vector field £ M with respect to which the 
"straight" differentiation operation $ is defined. 



3.3 * "Straight" differentiation of covariant fields 

The "straight" displacement scheme set up in the prececding subsection will 
evidently determine corresponding transformations of any continuous fields that 
may be given on the background space. The simplest case is that of an ordinary 
scalar field <P say, for which the (affinely parametrised) displacement x^ i— > x?^ 
will determine a corresponding pullback mapping <P i— > <P{ c -j with ^{ £ }{a;} = 
Evaluating this using an expansion of the form 

*{e) = # + *,,.(*&> - + l^u(^e } ~ *"Me] ~ + ^{c 3 } , (139) 



it can be seen from (135) and (137) that up to second order the required result 



will be given by an expansion of the standard form 

<P {e} = <P + e$$ + ^$ 2 <£ + 0{e 3 } (140) 
in which the relevant "straight" derivative coefficients are given by 

= $ 2 # = £T%V„<£ . (141) 

It is to be noticed that the formulae obtained in this way are of a strictly tenso- 



rial nature, unlike those in (137), whose non tcnsorially coordinate dependent 
nature is due to that of the coordinates to which they apply. It is also to be 
noticed that the formula for the second order "straight" derivative $ 2 <£ does not 
involve derivation of the displacement vector field £^ , unlike the formula for the 
corresponding second order Lie derivative, which is given by the expression £ 2 ^ 
= ^V M (fV^) = $ 2 <Z> + {pVpOVvf. 

After having thus dealt with the case of a scalar, the next simplest case to be 
considered is that of what Cartan would call a one- form, meaning just a covector, 

say, such as the electromagnetic gauge potential, for which the displacement 
x^ i — ► x { e \ wm determine a corresponding pullback mapping i— > A^^ with 

A{ e }fj,{%} = Av{ x {<i}} x { t } M - I n order to work this out, we need to evaluate the 
relevant tensor transformation matrix i,", whose components are the partial 
derivatives of the displaced coordinates with respect to the initial coordinates, 
so that they will be given by an expansion that is obtainable from (|134|) by 



partial differentiation of the coefficients (137) in the form 

_ e 2 



xA x = <*£ +e£\- ^(r^ P ,Are p + 2T^ p ^e\a) + 0(e 3 ) , (142) 
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in which it is of course to be understood that the fields and their partial coordi- 
nate derivatives indicated by the use of a comma, are all to be evaluated as the 
corresponding functions of the initial (undisplaced) coordinates x M . The cov- 
ectorial analogue of the scalar "straight" transportation formula (139) is thus 
obtainable in the standard form 

A^y^A^ + eU^ + ^A^ + Oie 3 } (143) 

with the relevant "straight" derivative coefficients given by 

$A P = rv„ + A v V»e , (144) 

$ 2 A M = rr(VpV^4 M - TV Va) + 2C(^ P )^A P , (145) 

where 7?. M „ A p is the Riemann tensor associated with the chosen connection, 
which is definable, using the (MTW) sign convention, by 

2V[ M V„]^ P = -n^ x p A x . (146) 

The cancelling out of all the non-tensorially coordinate dependent contributions 
so as to give a strictly tensorial result is of course not an accident but an 
automatic consequence of the covariant character of the definition of the straight 
diffcrcntation procedure. 

The next case that is of interest in the context of the kind of variational 
analysis that has been discussed in proceeding sections, is that of a two-index 
covariant tensor, B pl ,, which, if it were antisymmetric, as in the example of 
the Kalb-Ramond potential that was considered in previous work, would be 
intcrpretable as a two-form in the sense of Cartan. Whether or not it has 
such an antisymmetry property, the effect on any two-index covariant tensor 
of the "straight" transportation x^ L i— ► x? e y will, as in the previous case, be to 
induce a corresponding pullback mapping B pv i— > £>{ e } MI , which will b e give n 
by B {t}pu {x\ = B p<J {x {e} }x { p e} fi x { 1 } l/ . Working this out, using ( pi|) , ( p7| ), 
and (142) as in the previous single index case, one can obtain the corresponding 
expansion in the standard form 

B {e }^ =B„ v + e %B pu + -§ 2 B p „ + 0{e 3 } (147) 
in which the required "straight" derivative coefficients work out as 

$£W = £ P V P + B p ^e + B^i" , (148) 

% 2 B py = (%V CT iV - n pp \B Xv - K vp \B pX ) (149) 

+ 2e , ((Vp£ A )V p B A „ + (V^ x )V p B pX ) + 2ZV(V Ai e , )V„r • (150) 

The most important application of the preceeding result is to the case in 
which B pv is to taken to be the spacetime metric itself, for which one thus 
obtains 

9{e}»» = 9»u + e %9pu + -^% 2 9pu + 0{e 3 } (151) 
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in which the relevant straight derivatives are given by the comparitively simple 
looking formulae 

$g^ = 2V ( ^) , % 2 g pv = 2(V ( ^)V„ ) £ P - ^^Tl^ , (152) 

of which the former is the same as the usual Lie differentiation formula while the 
latter agrees with a result that was obtained in a different context by Boisseau 
and Letelier @. 

3.4 *Perturbed Dirac Goto Nambu action 

The calculus developed in the preceding sections is intended for application 
to cases in which one is concerned with a perturbed configuration that is ob- 
tainable from a given unpertubed perturbation by a small deviation, of order 
d say, via a continuously differentiable (nomotopic) variation, which may be 
taken to be specified by a homotopy parameter A say ranging from in the 
unperturbed configuration to 1 in the particular perturbed configuration under 
consideration. Using an overhead prime to denote differentiation with respect 
to A the "straight" transport vector £ M characterising the relative displacement 
of the final perturbed configuration (for which A = 1) will be expressible by an 
expansion of the form 

e = & + + 0{6 3 } . (153) 

Specifying the differential operator d with respect to the same parameter A, so 
that for example in the case of the action X we obtain the identification bX = X, 
the operator A giving the deviation between the original and perturbed value 
of any additive quantity will be expressible by a corresponding expansion 

A = S + \8 2 + <D{5 3 } . (154) 

The reason why we are interested in going to second order in the present 
notes is that - according to a very general principle first drawn to my attention 
many years ago by Taub [42]] , ]4^| - if a system of dynamical equations is obtain- 
able by a variational principle from an action X (meaning that the dynamical 
equations are the condition for the first order locally perturbed action 5X to 
vanish) then the corresponding first order perturbed dynamical equations will 
be similarly obtainable from the second order perturbed action S 2 X. 



For a free brane brane action of the kind governed by (131), the first order 
perturbed action with respect to a Lagrangian reference system (meaning one 
that is comoving with the brane worldshect) is directly obtainable (just by 
setting the external field variables and £> M „ set to zero) from (|88j) in the 
form 

51=^JdST^5 gilu , (155) 

and so the corresponding second order perturbed action will be obtainable, using 
(|9l"l), in the form 

5 2 X= l -JdS ({fnt9^+T^(\ir{5g pa )5g^ +S 2 g^)) . (156) 
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In the case to which our attention will be restricted here, the background 
metric will be taken to be held fixed in the sense that its Eulerian variation 
vanishes, Ag^ = 0. This means that according to the general displacement rule 



(138) the Lagrangian (comoving) metric variation in which we are interested will 
be given by 

A = ff$ + |(e$) 2 + <D{6 3 } . (157) 
The concise abbreviation $ for the straight diffcrcntation operator involved in 



(138) has been replaced here by the more explicit notation £$ to indicate its 



dependence on the displacement vector whose expansion is given by (153). It 
can thus be seen that the first and second order Lagrangian variations of the 
metric will be given by 

Sg MV = ('h^ > $ 2 g^ = (e$) V + ■ (158) 



When the relevant unperturbed field equation, namely (131), arc satisfied, 
the first order action variation due to any perturbation with compact support 
will of course vanish by the variation principle, 5 2 I = 0. In these circumstances 
the contribution from the second order displacement contribution £^ will auto- 
matically cancel out, so that for the Dirac-Goto-Nambu model determined by 
a surface Lagrangian of the trivial constant form 

< d >£ = -m d (159) 



t he s econd order action variation can be obtained from (156), using the formula 



(152) in the form 

s 2 z= j\J dsn(v(v^)v , r + z(y\p£ p )%]t u - ^;J p &) ■ (m 

with 

("ff^-m^Y*. (161) 

where to is a constant having the dimension of mass (which would be of the 
order of magnitude of the relevant Higgs mass scale in the case of a vacuum 
defect arising from the spontaneous symmetry breaking mechanism of the kind 
most commonly considered ||) and d is the dimension of the worldsheet (i.e. 
d=l for a simple point particle, d=2 for a string, and so on). 



Substituting the expression (161) for the surface stress energy momentum 
tensor into (131), the unperturbed Dirac Gotu Nambu equation of motion is 
obtained in a form that is given - independently of the dimension d indicated by 
the prefix ( d ) which may therefore be dropped - by the well known harmonicity 
condition that is expressible as the vanishing of the curvature vector, 

K» = 0. (162) 

In this rather degenerate special case, since there are no internal fields the 
linearised perturbation equation (|l32j) will of course completely determine the 
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evolution of the displacements by itself. By substituting (161) in (132) the 
corresponding perturbation equation is obtained in the form 

A (V^ A - rT^VX) = 2iV Vf . (163) 

It is straightforeward to check that this does indeed agree - subject to ( |162[ ) 
- with what is obtained [Ell by applying the variation principle to the second 



order action ( |160p . As well as the compact tensorial version (163) given here, the 
literature includes other equivalent but formally more complicated expressions 
[HI i Bjl (involving reference to int erna l coordinates or surface adapted frames 
of the kind discussed in Subection [bl]) that generalise earlier work which was 
restricted to a flat or De Sitter background [Et], |48| , fjgfl or specialised to the 
hypersurface supported ( "wall" or membrane) special case |50[| . 

3.5 *Higher order geodynamic models 

For most practical physical purposes the most useful generalisations of the Dirac- 



Goto-Nambu models governed by (162) and (163) are those of the very general 



category governed by (131) and (132) which allow for internal fields such as the 
currents that can be used, not only to represent the Witten type superconduc- 
tivity effect in cosmic strings, but also also to represent the effect of ordinary 
elasticity in terrestrial applications, such as the strings of musical instruments 
which were already a subject of scientific investigation, albeit at an empirical 
rather than theoretical level, in the time of Pythagoras. However before pro- 
ceeding to the discussion of such internal field effects in the following sections, it 
is of interest to consider how the simple Dirac-Goto-Nambu model can be gen- 
eralised in a way that goes beyond the "pure" brane description characterised 



in the free case by (131) and in the presence of an external force by (122). 

The distinguishing property of a "pure" brane model of the kind considered 
in Section ^| is the condition that the action depends only on the undifferentiated 
background fields <? M „, .4 p , B^ v , but not on their gradients. One of the most 
familiar kinds of example in which this condition fails to hold is that of an 
electromagnetically polarised medium, whose action |5l| depends not just on 
the gauge field but also directly on the associated field F^ u itself. For 
lack of time and space we shall not consider such electromagnetic effects here, 
but will consider only the simplest category of "geodynamic" brane models, 
meaning those in which, as in the "pure" brane models of the Dirac-Goto-Nambu 
category, the action depends only on the imbedding geometry of the worldshcet 
and not on any other external internal fields. The simplest such extension of 
the "pure" Dirac-Goto-Nambu model (whose action is proportional just to the 
surface measure which depends only on g^ but not its derivatives) is based on a 
Lagrangian consisting not only of a constant term but also of terms proportional 
to the two independent scalars that can be constructed as quadratic functions of 
the first derivatives of the metric, namely K^K^ and K ^„ p K^ 1 ,p ' . The inclusion 
of such "stiffness" terms has been suggested by Polyakov |32| and others, one 
of the main reasons being allowance for the deviations from the "pure" Dirac- 
Goto-Nambu description of cosmological string J53[, pM, fl55[] or domain wall 
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[ p6| , (57fl , Q , |39| , Q defects that one would expect to arise if the curvature 
becomes too strong. The kind of Lagrangian constructed in this way, namely 

< d >£ = -m d + b < d >A' p WK p - c ^K^ p WR^ P , (164) 

where m, b, and c are constants, has recently been the subject of investigations 
by a number of authors |6l| , |62) , ||| , Q , |65| , 66 . It is convenient to use the 
abbreviation 

(d)^P = b iA)rj t*» wk p _ c wk h*p j ( 165 ) 



which enables the Lagrangian (|16j) to be expressed in condensed form as 

< d >£ = -m d + WE^p WK^p . (166) 

The variation needed for evaluating the change in such an "impure" brane action 
will thereby be obtainable from the formulae above in the corresponding form 

|| 7 |r 1/2 £(||7l| 1/2 £) = {C+ l -C^5g, v (167) 

(again dropping the explicit reference to the brane dimension d) with 

SC = (£ xp »K A ; - 2 E^ p K x " p )5g^ + (2 - Z^ x )V x Sg^ . (168) 

Although it is still possible to construct a formally symmetric stress momentum 
energy density tensor of th e distributional type, the presence of the gradient 
term on the right of (168) will make it rather pathological, with not just a 
contribution proportional to a Dirac distribution as in the "pure" brane case 
described by (114) but also with a contribution proportional to the even more 
highly singular gradient of a Dirac distribution Q . In order to be able to con- 
tinue working in terms of strictly regular surface supported field, it is necessary 
[ p5| to deal with the gradient dependence of the action in such an "impure" 
brane model by having recourse to the use of a total stress momentum energy 
tensor T^ v of the no longer no longer symmetric canonical type that can be 
read out from the variation formula 

!l7ir 1/2 <S(||7l| 1/2 £) = T\V^ V + 2S l TTl p / v e - %(2Ef"V p C) (169) 



that is obtained after substitution of ( |126|) in (167) and ( |168|) . The regular but 
non-symmetric canonical surface stress momentum energy density tensor T p v 
obtained in this way []65| is given by 

T»„=f» v + t» v , (170) 

where T pv is symmetric and purely tangential to the worldsheet (as a regular 
"geometric" stress momentum energy density tensor would be) with a merely 
algebraic dependence on the second fundamental tensor, being given by 

f p v = ff = r\n\ = Zrf v - 2 S\ p K» p , (171) 
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while the remainder, which is of higher differential order, is expressible in terms 
of the third fundamental tensor (^) as 

t"„ = = 2 c S \» v - 2 b ~\\ , (172) 

which can conveniently be rewritten with the higher derivative contributions 
regrouped in the form 

K = 2(c - b)5» v +crf p r 1 ^n p aTX ±i , (173) 

where 

S% = ~\\ = ^ . (174) 

It is to be noticed that the total canonical surface stress momentum energy 
density tensor obtained in this way is is still automatically tangential to the 
worldshcct on its first (though no longer on its second) index, i.e. 

4J"„ = 0, (175) 

and that the higher derivative contribution proportional to the trace of the 
third fundamental tensor will drop out if the coefficients have the same value, 
c = b. 

The application of the variation principle to the effect that the surface in- 
tegral of the variation ( |169[ ) should vanish for any displacement £ M within a 
bounded neighbourhood can be seen to lead (via an application of Green's the- 
orem as in the Section ^| to dynamical equations |65[] of the form 

%T» v =2Zr°K p / v . (176) 

As in the pure Dirac-Gotu-Nambu case discussed above, the foregoing system 
of equations is partially redundant: although it involves n distinct spacetime 
vectorial equations, only n-d of them are dynamically independent, namely 
those projected orthogonally to the d-dimcnsional worldshcct. The others are 
merely Noether identities which follow independently of the variation principle 
from the fact that a displacement that is purely tangential to the worldshcct 
merely maps it onto itself and thus cannot affect the action, as can be verified 
directly using the generalised Codazzi identity (p5|). 

The lack of symmetry of T^ v means that the construction of the correspond- 
ing momentum current vector, V 11 say, associated with a generic background 
spacetime Killing vector field will not be quite as simple as in the "pure" 



brane case for which an expression of the simple form (119) suffices. However 



using the well known fact that the Killing equation (118) entails the intcgrability 
condition 

V^P = Tl p upX k x , (177) 
together with the observation that the antisymmetric part of the canonical stress 



momentum energy density tensor (171) is given according to (172) just by 



r [H = _2V A £ A[M,yl , (178) 
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it can be seen that, for any solution of (118), the ansatz p 



V p = T* v k v + 2S^ v p V v k p (179) 
provides a surface current, , which satisfies the tangentiality condition 



J" 

and for which the strict surface conservation law 



A -P M = 0, (180) 



V^=0, (181) 



will hold whenever the equation of motion (176) is satisfied. It is to be remarked 
however that for a Killing vector of the irrotational kind for which V^fc^ vanishes 
altogether the second term (interpretable as a surface spin density contribution) 
in ( |179| ) will not contribute, i.e. an expression of the simpler form ( |119| ) will 
suffice. This applies in particular to the case of an ordinary translation gen- 
erator in flat space, for which the corresponding conserved surface current will 
represent ordinary energy or linear momentum, whereas in the case of angular 



momentum the extra (spin density) term in (179) is indispensible 



The proceeding formulae all include allowance for arbitrary background cur- 



vature, but, to obtain the analogue of the non redundant version (131) of the 
equations of motion in a reasonably simple form, the restriction that the back- 
ground spacetime be flat, i.e. 1Z p up<J = 0, will now be imposed. This enables 
the required system of dynamical equations to be expressed |65|| in the form 

f^K^p = 2(b - c)r p % (J? V V" A- a ) , (182) 

with the higher derivative terms grouped on the right hand side, which vanishes 
if b = c. 

In the particular case [6C[| of a membrane, meaning a brane supported by a 
hypersurface of dimension d=n-l, the second fundamental tensor and its trace 
will be given in terms of the unit normal X p (which in this case will be unique 
up to a choice of sign) by A' M „ p = K^ v \ p and K p = K\ p with K = K v v where 
is the second fundamental form (whose sign depends on that orientation 
chosen for the normal). In an ordinary 4-dimcnsional spacetime background, 
this membrane case corresponds to d=3, for which the symmetric tangential 
part of the surface stress momentum energy density tensor will be therefore be 



obtainable from (171) in the form 



WW" = -(m 3 -b< 3 lA 2 -K (3) AV 3) A°p) ( 3 '7^-2b< 3 >A< 3 »A' il/ + 2c< 3 )A A ' i < 3 )A A I/ . 

(183) 

In the case of a string with d=2 one can use the fact that the trace free confor- 
mation tensor @ will satisfy 2 ( 2 >C A pp < 2 >C Xl f = WC K \ P mC KXp Mr]^ to obtain 
a corresponding formula (which holds regardless of the background spacetime 
dimension n) given |65| ] by 

wf» v = -m 2 V>rf v + 2(c - b) ( 2 >A' p . (184) 
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4 Strings and other spacially isotropic brane mod- 
els 



4.1 The general category of "perfect" brane models 

It is reasonable to postulate that a "weak" energy condition of the kind formu- 
lated and justified by Hawking and Ellis [|l5| should hold for any pure p-brane 
model as a condition for physical realism as a macroscopic description of a 
(p+l)-surface supported physical system at a classical level, meaning that the 
model's surface stress momentum energy density tensor (as introduced in Sec- 
tion ||) should be such that the contraction T^P^fi^ is non negative for any 
vector /3 M that is timelike. Furthermore the causality condition to the effect 
that there should be no timelike characteristic covector (i.e. no supcrluminal 



propagation) can be seen from (133) to entail the further requirement (going 
marginally beyond the "weak" condition of Hawking and Ellis) that T^P^P^ 
should be strictly positive if P" is timelike. This leads to the formulation of 
what may be called the "minimal" energy condition for a pure p-brane which is 
expressible as 

& t p lt <0 T^P„p u >0. (185) 

A (pure) p-brane model will consequently be characterised by a well defined 
surface energy density, U say, that is specifiable by an eigenvalue equation of 
the form 

T^vfi 1 ' = -Up" , (186) 

where the corresponding eigenvector /3 M is distinguished by the requirement that 
it be tangential and non-spacclike: 

V=0, /3"/3„<0. (187) 



It is apparent that the "minimal" energy condition (185) requires that the eigen- 



value U should be strictly positive unless P" j s null in which case it may vanish: 

U>0, U = => /^/3 M = 0. (188) 

In the Dirac-Goto-Nambu model that is most familiar to present day cos- 
mologists, the eigenvector /3 M is indeterminate and the energy density U is the 
same (in relativistic units such as are used here, with the speed of light set to 
unity) as the corresponding (surface) tension. However in general it is essential 
to distinguish the concept of energy density U from the concept of the tension 
(as used in physics since the formulation of Hooke's law at the time of Newton) 
from which term "tensor" is derived. The (surface) tension scalar, for which we 
shall use the traditional symbol T, is defincable gcncrically in a manner con- 
sistent with traditional usage, for a "pure" p-brane (i.e. a (p+l)-dimensionally 
supported system) of the kind considered here by decomposing the trace of the 
surface stress momentum energy density tensor in the form 

= -U - pT . (189) 
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Apart from the degenerate Dirac-Goto-Nambu case for which U and T are 
actually equal, the simplest possibility is that of a perfect p-brane || , meaning 
one whose surface stress momentum energy density tensor is spatially (thus p- 
dimcnsionally) isotropic so that it is will be expressible, for a suitable choice of 
the normalisation of the eigenvector /3 M in (|186|) by 



T p,v = pup _ Trj nv : ^gg) 

where the required normalisation is given by 

U — T = < . (191) 



This category includes the case of the Dirac-Goto-Nambu model (161), which 
is obtained (with T = mP +1 ) by normalising the (in this case indeterminate) 
eigenvector to zero, i.e setting /? M = 0. A more mundane example is provided 
by the familiar "improper" case p=n-l (where n is the background dimension) 
of an ordinary (rclativistic) perfect fluid with pressure P = —T. Although the 
tension T is negative in the ordinary fluid case, it must be positive, as a con- 
dition for stability, for "proper" p-branes of lower dimension, p<n-l. Whereas 
the membrane models (with p=2 for n=4) that are appropriate for the descrip- 



tion of nautical sails will not in general have the isotropic form ( 19C ) , a familiar 
everyday example of a membrane that will automatically be "perfect" in this 
sense is provided by the case of an ordinary soap bubble whose boundary hy- 
persurface will be characterised by < T <C U. The kind of application on 
which we shall concentrate in the following sections is that of a string, as given 



by p=l, for which the "perfection" property (190) will always hold: isotropy 
cannot fail to apply in this 1-branc case because only a single space dimension 
is involved. 

Before specialising to the (automatically perfect) string case, it is to be re- 
marked that except when the /3 M is indeterminate (as in the Dirac-Goto-Nambu 
case), or null (as can occur in special subspaces where the current becomes 
null in a string model of the kind ||(|, |37]], (4lJ appropriate for describing cos- 
mic vortex defects of the "superconducting" kind proposed by Witten ||) this 
eigenvector will determine a corresponding unit eigenvector, 

u" = (C/-T)~ V V, u^u M = -l, (192) 



specifying a naturally preferred rest frame, in terms of which (19C) will be 
rewriteable in the standard form 

T"" = (U -T)«V -Trf" . (193) 

With respect to the frame so defined, it can be seen that the velocity c B say 
of propagation of extrinsic perturbations (what are commonly referred to as 
"wiggles" ) of the perfect branc worldshcct will be given, according to the char- 



acteristic equation (133), by 

T 



c * 2 = Tj- 
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By substituting (192) in the (125), and using the expression (18) for the accel- 



eration of the preferred unit vector (192), the extrisic equation governing the 



free motion of the support surface of any (proper) perfect brane is reducible to 



the standard form M, |41 



cjK» = (1 - c E 2 )^ . (195) 

In the Dirac-Goto-Nambu case one has U — T and hence c E 2 = 1, (i.e. the 
"wiggle" speed is that of light) so that the right hand side of (|l95| ) will drop 
out, leaving the harmonicity condition that is expressed by the vanishing of the 
curvature vector given by (|2l|). At the opposite extreme from this relativistic 
limit case, one has applications to such mundane examples as that of an ordinary 
soap bubble membrane or violin string for which one has c E 2 <C 1. 

4.2 The special case of "barotropic" brane models 

Except in the Dirac-Goto-Nambu limit case, the extrinsic equation of motion 



(196) will need to be supplemented by the dynamic equations governing the 
evolution of the internal fields, and in particular of T, U, and u^, on the brane 
surface. The simplest non trivial possibility is the case of a perfect brane model 
that is "barotropic" , meaning that its tension T (or equivalently, in the case of a 
fluid, its pressure P) is a function only of the energy density U in the preferred 
rest frame. In this barotropic or "perfectly elastic" case, the complete set of 
dynamical equations governing the internal evolution of the perfect brane will be 



provided just by the tangential projection of the force balance law ( |110| ) which, 
in the free case to which the discussion of the present section is restricted, has 
the simple form 



V X ^T^ = 0. (196) 

In the barotropic case (which includes the kind of cosmic string models ]36| , 
[ p7[ , Q appropriate for describing the Witten type superconducting vaccuum 
vortices Q whose investigation provided the original motivation for developing 
the kind of analysis presented here) experience of the "improper" case of an 
(n-l)-brane, i.e. the familiar example of an ordinary barotropic perfect fluid, 
suggests the convenience of introducing an idealised particle number density, v 
say, and a corresponding (relativistic) chemical 'potential, or effective mass per 
idealised particle, \i say, that are specified modulo an arbitrary normalisation 
factor by the equation of state according to a prescription of the form 

which fixes them modulo a pair of constants of integration which are not al- 
lowed to remain independent but are related by the imposition of the restraint 
condition 

fj,u = U-T. (198) 
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The quantities introduced in this way can be used for defining a surface particle 
current density vector v p and a corresponding dynamically conjugate energy - 
momentum (per particle) covector fx p that are given by 

v p = vu p , jip = (j,Up , (199) 



in terms of which the generic expression (193) of the surface stress momentum 



energy density tensor of a perfect brane can be rewritten in the form 

T\ = z^/v - Trf a , (200) 
whose divergence can be seen to be given identically by 

%TP a = ^%v p + 2z/%/V] - TK a . (201) 
In the free case for which this divergence vanishes, projection orthogonal to the 



worldshcct gives back the extrinsic dynamical equation (195), while by contrac 



tion with v p one obtains the simple surface current conservation law 

%v p = . (202) 
The remainder of the system of equations of motion for the free perfect brane can 



hence be obtained from (196), as a surface generalised version of the standard 



perfect fluid momentum transport law J67[ , in the form 

"V^/ip] = , (203) 

which is such as to entail a corresponding Kelvin type law of conservation (under 
Lie transport by the flow) for the relevant surface vorticity tensor if^^J^Hp] 
so that if it vanishes initially it will remain zero throughout, which is the inte- 
grability condition for the existence of a surface potential function, ip say, such 
that [ip = V p i/?. 

Appart from the extrinsic perturbations of the world sheet location itself, 
which propagate with the "wiggle" speed c E (relative to the frame deterined 
by u p ) as already discused, the only other kind of perturbation mode that can 
occur in a barytropic string are longitudinal modes specified by the variation of 
U or equivalcntly of T within the world sheet. As in the "improper" special case 
of an ordinary perfect fluid, for which they are interpretable just as ordinary 
sound waves, such longitudinal "woggle" or "jiggle" perturbations can easily be 



seen from (U00h and (201) to have a relative propagation velocity c L say that is 
given D, g| by 



2 v dfi dT 



- Hdv dU- (2 ° 4) 

In order for a barotropic model to be well behaved, c B 2 and c L must of course 
both be positive, in order for the velocities to be real (as a condition for local 
stability), and they must also both be less than unity in order to avoid supcrlu- 
minal propagation (as a condition for local causality). An important qualitative 
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question that arises for any proper barotropic p-brane, i.e. one with 1< p<n-l 
(excluding the trivial case p=0 of a point particle, and at the other extreme the 
"improper" case p=n-l of a fluid) is whether the extrinsic propagation speed 
is subsonic, meaning that "wiggles" go slower than "jiggles" , i.c c E < c L , tran- 
sonic, meaning that the speeds are the same, c E = c L or supersonic, meaning 
that "wiggles" go faster than "jiggles", c E > c L . Much of the earliest work p8| |, 
H> ®> 0. [@, @ on the Witten type superconducting vaccuum vor- 
tex phenomena was implicitly based on use of a string model of subsonic type, 
but a more careful analysis of the internal structure of the Witten vortex by 
Peter [frs) , [f76| , |nj has since shown that models of supersonic type are more 
appropriate. On the other hand such everyday examples as that of an ordinary 
violin string are subsonic, a condition that is sufficient |7^] though not neces- 
sary |79| , |8(J for stability of the corresponding circular ccntrifugally supported 
"vorton" type loop configurations |8l| , Q , Q , |84|] that will be discussed later 
on. 

A particularly important special case is that of a barytropic p-brane of per- 
manently transonic type, meaning one characterised by an equation of state of 
the non dispersive constant product form 

UT = to 2p+2 , (205) 

where m is a constant having the dimensions of a mass, for which application of 



the formulae (194) and (202) gives c E 2 = c L 2 not just for some critical transition 
value of U or T (as can occur for other kinds of equation of state) but for all 
values. In the string case, p=l, it can easily be demonstrated J5| that this 
non dispersive permanently transonic model represents the outcome of a (rather 
artificial) dimensional reduction process first suggested by Nielsen Q , Q , ]87j , 
[ p8| . As will be described below, this model can be shown J89| to be governed by 
equations of motion that are explicitly integrable in a flat empty background, 
and it can be used (89) , |)0j to provide what (contrary to a misleading claim 
that has been made) can in fact (93) be a highly accurate description of 
the macroscopic effect of "wiggles" in an underlying cosmic string of the simple 
Goto-Nambu type. 



5 Essentials of elastic string dynamics 

5.1 *Bicharacteristic formulation of extrinsic equations 

Between the hypcrsurfacc supported case of a membrane and the curve sup- 
ported case of a point particle, the only intermediate kind of brane that can 
exist in 4-dimensions is that of 1-brane, i.e. a string model, which (for any 
background dimension n) will have a first fundamental tensor that is express- 
ible as the square of the antisymmetric unit surface element tensor E^ v = —£ u>1 



given by (|75j): 



rf v = 8%£» v . (206) 
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In the case of string (to which the remainder of this article will be devoted) the 
symmetric surface stress momentum energy tensor T^ v that is well defined for 
any "pure" (i.e. unpolarised) model will not only be expressible in the generic 



form (190) that expresses the (in the string case trivial) property of spacial 
anisotropy with respect to the preferred timclikc or null eigenvector /3 M , but it 
will also be expressible generically in bicharacteristic form as 

= (3^f3 u ] , (207) 

in terms of a pair of timelike or null tangent vectors (3*? that can be seen to 
be extrinsic bicharacteric vectors, meaning that they lie respectively along the 
two ( "right moving" and "left moving" ) directions of propagation of extrinsic 
perturbations, and are thus respectively ort hogonal to the covectors x M given 



by the extrinsic characteristic equation (133) whose solutions can be seen from 
( |207|) to be given by 0£xn = or /3 M x M = 0. The expression (|207D docs 
not completely determine the bicharacteristic vectors /3£, but leaves open the 
possibility of a reciprocal multiplicative rescaling of their magnitudes, subject 
to the condition that their scalar product remains invariant, its value being 



obtainable by comparison of (EOT]) with (190) as 



P +Ii p M = T v v = - (U + T) . (208) 

The geometric mean of the magnitudes of the bicharacteristic vectors remains 
similarly invariant, with value given by the magnitude of the preferred eigen- 
vector Z?^ as given by ( |191[ ), i.e. 

(/V/?/) {P-uP-) = {P»P? = (U - T) 2 . (209) 
To obtain the equations of motion of the string we need to use the force 



balance equation (112), which just takes the form 

%T" V = T , (210) 

for a string that is isolated, where / M (which simply vanishes in the free case) 
is the force exerted by any background fields, such as the electromagnetic field 
and a mbient (Higgs or other) fluid that are allowed for in the explicit expression 



( [1131) . (The only difference if the string were not isolated, i.e. if it belonged 



to the boundary of one or more attached membranes, is that an additional 



contact force contribution of the form (|112| ) would also be needed as well as 
the background contribution / M .) Using the Wcingarten intcgrability condition 
( p"5| ) which is equivalent in the string case to the the projected Lie commutativity 
condition 

SE v {fi/V ll p_ v - PJ*W + V ) = , (211) 



it can be seen from ( p07| ) that the extrinsic equations of motion governing the 
evolution of the string worldsheet will therefore &J be expressible in character- 
istic form as 

4^/3/V^/ = J*,/" , (212) 
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for either choice of the sign ±, one version being obtainable from the other and 
vice versa by (211). 

Except where the string state is locally degenerate in the sense that one of the 
bicharacteristic vectors (3 J 1 is null, which by ( plOj ) will only occur where T = U 
(as is everywhere the case for the Goto-Nambu model, for which both bicharac- 
teristic vectors are null everywhere) it will be possible to fix their normalisation 
by requiring that they have the same magnitude, /3 ±M /? ± M = —{T — U) and 
hence to define a corresponding pair of unit bicharacteristic vectors 

u ± *= (U-Ty 1/2 (3 ± » , u ± ^ = -l. (213) 

5.2 Preferred orthonormal diad for nondegenerate state 

In the generic case, T < U, for which the eigenvectors of the stress-energy 
tensor are non null, the corresponding unit vectors will constitute a preferred 
diad consisting of the unit timelike eigenvector u M introduced in (192), together 
with an orthogonal spacelike unit eigenvalue v^. These mutually orthogonal 
unit vectors will be related to the corresponding unit bicharacteristic vectors 
u M and u_ M (as given by (|213|) above) by the expressions 



u" = v/l 2 Ce2 -(u + ^ + u_») , t,M = V^^_( u+ /*_„_/*). (214) 

In terms of this preferred diad, the antisymmetric unit surface element E^ v 
and the (first) fundamental tensor rf" v of the string worldshcet will be given by 

= u^v v - v"u v , if = -u^u" + v^v" , (215) 



and the standard form (193) for the surface stress momentum energy density 
tensor will reduce to the simple form 

T" v = Uu^u" - r// . (216) 



Proceeding directly from this last form, t he ex trinsic dynamical equation ( |211| ), 
i.e. the surface orthogonal projection of (|210| ), can be rewritten as 

±%(Uu" -Tv'")=±^r , (217) 

using the notation 

u» = u v W v u v , = v v V u v v . (218) 

in terms of which the worldsheet curvature vector is expressible as 

= A? v (v ,v - ii v ) . (219) 
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5.3 Duality in elastic string models 

The formulae in Subsection [T^ illustrate a special feature distinguishing string 
models from point particle models on one hand and from higher dimensional 
brane models on the other, namely the dual symmetry |3(| that exists at a formal 
level between the timelike eigenvector and the associated eigenvalue U on 
one hand, and the the dynamically conjugate quantities that are the spacelike 
eigenvector and v M and the associated eigenvalue T . This feature is particularly 
informative in the barotropic or "perfectly clastic" case discussed at the end of 
the Section |], i.e. in the case when the internal dynamics is controlled just by 
an equation of state giving the energy density U as a function of the tension T 
only, since in there will be a corresponding duality relation between the number 



density v and the associated chemical potential \i as introduced in (197), which 



is equivalent to the mutually dual pair of differential defining relations 

dU = [xdv , dT = -vdn . (220) 

Defining the surface duals of the current vector v p and momentum covector jj, p 



given by (199) as 

*v p = 8 pa v° = vv p , *// = £ = fiv" , (221) 



one can use (201) to obtain the force balance equation for a barotropic string 



model in the self dual form 

M P V CT ^ + ★i/pVa**" + J?„ (Uu a - Tv" 7 ) = J p , (222) 
whose orthogonally projected part gives back the self dual extrinsic force balance 



equation (217), while its tangcntially projected part gives a mutually dual pair 



of surface current source equations of the form 

(U - T)%v° = -v~j p , (17 - T)%ip° = *fj, p J P ■ (223) 

It can be seen from this that, like the surface number current density v a , the 
dual current */i cr will also be conserved in a free barotropic string, i.e. when the 
external force f p is absent. 

It is useful pq| to introduce a dimensionless state function called the "char- 
acteristic potential" d say that is constructed in such a way as to be self dual 
according to the differential relation 



d# = J^. (224) 

V l lv 

This function is convenient for the purpose of writing the internal equations of 
motion obtained by tangential projection of ( |210|) in a characteristic form that 
will be the analogue of the characteristic version (|212|) of the extrinsic equations 



of motion obtained by the orthogonal projection of (210). The longitudinal 



bicharacteristic unit vectors IJ 1 say (the analogues for the internal sound type 
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waves of the extrinsic bicharacteristic unit vectors uj 1 introduced above) are 
dcfincable by 



f/ = (l-c L 2 ) V V±ClO , W = -l. ( 22 5) 



where c L is the longitudinal characteristic speed as given by (204). In terms 
of these, the tangentially projected part of the force balance equation ( |210| ) is 
convertible |36[ ] to the form of a pair of scalar equations having the characteristic 
form 

(U - T)£/ (V p ti T v v %u v ) = ±£^ £/ J" , (226) 
from which it is directly app arent that, as asserted in the Section |J, the state 



function c L , as given by (|204| ), is indeed the characteristic speed of longitudinal 
"woggle" propagation relative to the preferred frame specified by the timelike 
eigenvector u M . When rewritten in terms of the unit bicharacteristic vectors uj 1 
that are the extrinsic analogues of £ ± M , the corresponding extrinsic equations 



(212) take the form 

(U - T)J>/V P = ±»J" , (227) 

(It is to be noted that the corresponding equation (32) at the end of the original 
discussion |36| of the characteristic forms of the equations of string motion 
contains a transcription error whereby a proportionality factor that should have 
been \f(JJ — T)T was replaced by its non-relativistic - i.e. low tension - limit, 
namely VUT.) 



5.4 *The Hookean prototype example of an elastic string 
model 

What essentially distinguishes different physical kinds of "perfectly clastic" , i.e. 
barotropic, string model from one another is the form of the equation of state, 
which depends on the microscopic internal details of the vacuum defect (or other 
structure) that the string model is supposed to represent. In a cosmological con- 
text, the most important special case, whi ch w ill be discussed in detail later on, 
is given by the constant product formula (|205| ), but for terrestrial applications 
the most useful kind of equation of state for general purposes is the one named 
after Robert Hooke, who is recogniseably the principal founder of string theory 
as a branch of physics in the modern sense of the word. (As for superstring 
theory, it remains to be seen whether it will become established as a branch of 
physics or just as a topic of mathematics.) Hookc's famous law to the effect 
that the tension, T is proportional to the extension, i.e. the change in \ jv in 
the notation used here, provides a very good approximation in the low tension 
limit that is exemplified by applications such as that of an ordinary violin string. 
This law is expressible within the present scheme as 

U{Y + T) = U a Y + \T 2 , (228) 
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where U is a constant interpretable as the rest frame energy in the fully relaxed 
(zero tension) limit, and Y is another constant interpretable as what is known 
(after another pioneer worker on the subject) as Young's modulus. 

5.5 General elastic string models 

A more general category of non-barotropic string models is of course required, 
not only in such terrestrial applications as (above ground) electric power trans- 
mission cables (the subject of my first research project as an undergraduate 
working working for commercial industry during a vacation) , but also for many 
cosmologically relevant cases that might be envisaged (such as that of a "warm 
superconducting string" j94|) in which several independent currents are present. 
Nevertheless most of the cosmological string dynamical studies that have been 
carried out so far have been restricted (not just for simplicity, but also as a very 
good approximation in wide range of circumstances) to models of the the "per- 
fectly elastic" barotropic type, and more specifically to the category (including 
the Hookean and Goto Nambu examples as opposite extreme special cases) of 
string models describable by a Lagrangian of the form 

Z = K + ~fA ll + \w^B IJlv , (229) 

in which the only independent field variable on the world sheet can be taken 
to be a stream function ijj, which enters only via the corresponding identically 
conserved number density current 

d 1 = , Vv = V„V , (230) 

whose squared magnitude, 

X = = -p» Pu = -f!^ , (231) 

is the only argument in a (gencrically non-linear) master function, 

A = A{ x } . (232) 

(It is to be noted that there was a transcription error whereby the factor 1 / 2 



in the final term of ( 229 ) was left out in the preceeding version of these notes 
[^.) In order to satisfy the vorticity flux conservation requirement (as derived 
in Section || from the condition of invariance with respect to gauge changes in 
the background Kalb-Ramond field B^) one must have 

= h£^ , (233) 

for some fixed (quantised) circulation value h say - a sort of helicity constant 
associated with the string. In the simplest "local" string models p(i| ], |37j this 
quantity h can simply be taken to vanish, while in the in the simplest "global" 
models p2] ], fl34| , |3f| it will be directly identifiable with Planck's constant, i.e. 
one will have h = 2nh. (In more general cases it might be an integral multiple of 
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this.) The electromagnetic current density is specified in terms of the gradient 
of the stream function in the form 



f = qcf , (234) 

where q is another fixed parameter which is interprctable as the electromagnetic 
charge coupling constant per idealised particle of the current, which means that 
its value will depend on the convention used to fix the normalisation of </' or 
equivalently of c M , a question that will be taken up in the discussion following 



equation (247) below. 

Previous studies of cosmic strings have mostly been restricted to cases involv- 
ing allowance for a Kalb-Ramond coupling (which is relevant for vortex defects 
of "global type" p^j , |33|j , |34{| , pq| , [Q in the absence of an internal current, 
or alternatively to cases (such as string models of the kind J|(|, || needed to 
describe local superconducting vortices arising from the Wittcn || ) mechanism 
involving allowance for an internal current in the absence of Kalb-Ramond cou- 
pling. However there is no particular difficulty in envisaging the presence of 
both kinds of coupling simultaneously as is done here. A new result obtained 
thereby (in |^) is that for stationary and other symmetric configurations, the 
effect of the coupling to the Kalb Ramond background is expressible as that of 
a fictitious extra electromagnetic field contribution. 

It is the specific form of the master function that determines the relevant 
equation of state - or if its range of validity is sufficiently extensive, the pair 
of equations of state, corresponding to the qualitatively different regimes dis- 
tinguished by the positivity or negativity of \ - by which the dynamics of the 
model is governed. The derivative of the master function provides a quantity 

K = 2^ , (235) 
dx 

in terms of which the surface stress momentum energy density tensor defined 



by (90) will be given by 

= Kp^p v + Arj^ . (236) 

This only depends on the internal field gradient and on the imbedding and the 
background geometry: it involves neither the background electromagnetic and 
Kalb Ramond fields nor their respective coupling constants q and h. The role 
of the latter is just to determine the b ackgr ound force density, which will be 



given, according to the general formula ( 113 ) , by 



f p = qF m c» + \n p ^» . (237) 



Unlike the force density, the stress momentum energy density tensor, as given by 
( [236|) , is formally invariant under a duality transformation 0], f36f whereby the 
variational momentum covector p p , the current c M , and the master function, 
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A itself, are interchanged with their duals, as denoted by a tilde, which are 
definable by 



c" = IC P P , P^JCc^ x = = -p% , (238) 

and by the relation 

XK, = A - A = - X IC , (239) 

with 

K = 2^=K~\ X=-IC 2 X- (240) 

dx 

In terms of these quantities the stress momentum energy density tensor can be 
rewritten in either of the equivalent mutually dual canonical forms 

c?p v + Atft = T"„ = &% + V„ • (241 ) 



Substitution of these respective expressions in (210) then gives the force balance 
equation in the equivalent mutually dual forms 

Pp %c p + 2c p % Pp] + KK p =J P = p P %c p + 2c^%p p] + AK P . (242) 

The independent internal equations of motion of the system are obtainable by 
contraction with the independent mutually orthogonal tangent vectors p p and 
p p , the latter giving what is just a kinematic identity, since, by construction 
according to (p37|) one automatically has 



P p f P = (243) 

and hence 

%c» = & £ ptl %p p = 0. (244) 

The only part of the internal equations of motion that is properly dynamical 
(from the point of view of the variation principle with respect to A) is thus just 
the part got by contraction with p p : since (237) gives 

P P 7 P = qX^ p F pp . (245) 

one obtains the equivalent alternative forms 

= g£ pp F pp & £ pp (2%p p] + qF w ) = . (246) 

It is to be noted that unlike the electromagnetic force, the Joukowski force 
(interpretable as a manifestation of the Magnus effect) arising from the Kalb- 
Ramond coupling in (237) will always act orthogonally to the world sheet and 
thus does not affect the internal dynamics of the string. 
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5.6 Standard normalisation of current and charge 

The preceeding equation ( |246| ) can be construed as the (Poincare type) intc- 
grability condition for the local existence on the string worldsheet of a scalar p 
in terms of which the dual momentum will simply be given as the tangentially 
projected gauge covariant derivative: 

p p = %D p p , D p = rf p D„ , Dp = V„ - eA h , (247) 

where tt is a normalisation constant and e is a new, implicitly more fundamental, 
charge coupling constant in terms of which the original charge per unit idealised 
particle of the current is given by 

g = jje. (248) 

The scalar tp can be used as the primary independent variable, instead of ip, in 
an alternative (dynamically conjugate) variational formulation 0, p6| ] using a 
Lagrangian constructed from A instead of A. 

By adjusting the scale of calibration used for the measurement of either the 
current d 1 or the field ip it would of course always be possible, for purposes of 
mathematical convenience to simply set the dimensionless scale factor tt to unity. 
However other choices may be more useful for purposes of physical interpreta- 
tion. In particular, for the kind of string model ^ that provides a macroscopic 
descrition of Witten type syperconducting vortices, the scaling of ip can most 
naturally be chosen in such a way as to allow this field to be interpreted as the 
phase of the complex scalar field characterising the underlying bosonic conden- 
sate, so that the latter will be expressible in the form <P = \<P\c lcp . On the other 
hand it has been found mathematically convenient in the same context to fix 
the normalisation of the current c M in a standard manner by requiring that the 



coefficient K, introduced in (235) should tend to unity in the null current limit, 
i.e. as x tends to zero. The condition for this standard current normalisation to 
be consistent with the natural normalisation of p is expressible as the require- 
ment that the ratio jj (interpretable as the number of fundamental particle units 
per idealised particle unit of the current c M ) should just be a square root of the 
model dependent constant k that was used in previous work Q, i.e. ft 2 = k . 
The corresponding fundamental particle current vector z* 1 say will thus be given 
by 

z" = tic'' , (249) 
so that if the particles are charged, with fundamental charge coupling constant 



e, the corresponding electromagnetic surface current vector (234) will be given 
simply by 

J" = ez>* . (250) 

If its range of definition is sufficiently extended, the master function A will 
determine not just one but a pair of distinct equations of state, one applying to 
the "magnetic" regime where the current d 1 is spacelike so that x is positive, 
and the other applying to the "electric" regime where the current cf is timelike 
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so that the \ is negative. In cither of these distinct regimes - though not on the 
critical intermediate "null state" locus where the current is lightlike so that \ 



vanishes - the mutually dual canonical forms (241) of the stress momentum en- 
ergy density tensor are replaceable by the equivalent more manifestly symmetric 
standard form 

T" v = -&c v + -cPV . (251) 
X X 



By comparison with the corresponding expression (216) it can thus be seen that 
the state functions introduced in the variational function that has just been 
specified will be phyiscally interpretable in terms of the corresponding energy 
density U and tension T, and also of the number density v and effective mass /i 
associated with the corresponding physically normalised particle flux (meaning 
that for which the charge per particle is e) will be given paramctrically in the 
"electric" regime where c M is timelike, by 

A = -17, A=-T, fx = ~v\ X = ttV, ^ = t~, (252) 

v 

with 

He' 5 = vu p , c p = %[iv p , %p p = vv p , p p = ftfiup . (253) 

(For the corresponding mathematically normalised flux, whose idealised parti- 
cles are to be considered as each containing jj of the fundamental particles - so 



that, in accordance with (248) the corresponding charge per particle is jje - the 
corresponding effective mass would be jj/i and the corresponding number density 
would be The analogous formulae for the "magnetic" regime where c p is 

spacelike (the first possibility to be considered in the early work on the Witten 
mechanism H), will be given by 



A = -T, A=-U, tx^ii\ X = -f" 2 , £ = tt 2 -, (254) 



/' 



with 



jjc p = \iv p , ~c p = $vu p , ((pp = fiUp , p p = %vv p . (255) 

Whether or not it is charged, and whatever the local (spacelike, timclike, or 
null) character of the current may be, the conservation law ( |236| ) implies that 
its integral round any closed string loop will give a corresponding conserved 
particle quantum number 

Z = ldx p £p V z v . (256) 



For such a loop there will also be a second independent topologically conserved 
phase winding quantum number, N say, defined by 

2ttN = <£d<p. (257) 



It can be seen from (247) that this number will be given by 

2n$N = g$ + I dx p P p , (258) 
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where $ is the magnetic flux integral, 



3> = jdx^A^. (259) 

It is to be noted that in the clcctromagnctically charged coupled case, i.e. when 
the relevant fundamental particle coupling constant e has a non zero value 
(whose approximate value will be given in fundamental units by l/\/137 or 
a multiple thereof) the loop will have a conserved total electric charge, Q say, 
given by 

Q = eZ = jdx^E^'f . (260) 

5.7 Analytic string model for Witten vortex 

In realistic conducting string models one can not expect the appropriate master 
function A{%} to be given exactly by an explicit analytic formula, but only by 
the output of a detailed numerical computation of the internal structure of the 
underlying vacuum defect Q , |Q , Q , (77) . Nevertheless it has recently been 
found H that a remarkably good approximated description of the results of such 
computations for Witten's original "superconducting" bosonic vortex model can 
be provided by analytic expressions of a fairly simple kind, involving only two 
model dependent mass parameters. The one that is presumed to be the larger is 
the usual Kibble mass to say, whose square determines the energy and tension 
in the null current limit, and whose value is expected to be of the order of the 
mass of the Higgs field responsible for the spontaneous symmetry breaking that 
produces the vortex defect. The other mass parameter, to* say, is expected 
to be of the order of the considerably smaller mass scale that is presumed 
to characterise the independent current carrier field in the relevant version of 
Witten's bosonic field model ||. In all the cases that have been investigated 
in detail |75[ j, |76[ ], ]77[ | it has been found that to obtain a reasonably accurate 
representation throughout the ("magnetic") regime for which x is negative, and 
for moderate positive values as well, it suffices Q to use the purely polynomial 
formula 

A =- ma+ !( i -2^)" 1 > (26l) 

which is applicable within the range 

6 2m* A m z + m* z 

which is limited below by the condition that c L 2 should be positive and limited 
above by the condition that c E 2 should be positive. For va lues near the upper 
end of this range (in the "electric" regime) the formula ( [26 1|) docs however, 
becomes inadequate, but to obtain a very satisfactory description throughout 
the positive range of x, it suffices to replace it by the still very simple though 
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non-polynomial formula 

A = -m 2 -^ln{l-^}, (263) 

for which the corresponding range of applicability is given by 

- 1 < -^r < 1 - e - 2m2 / m * 2 . (264) 

Since the underlying field theoretical model proposed by Witten || provides 
what can at best be only a highly simplified description of the mechanisms that 
would be involved in any real vacuum vortex, there is no point in investing too 
much effort in describing it with unduly high precision. To provide just the 
moderate degree of accuracy that is physically reasonable, the single formula 
( |263|) is all that is needed to cover the entire admissible range: it is more accurate 
than it need be for this purpose in the positive range of x; it agrees perfectly 



with the preceeding formula (261) in the null limit where x tends to zero; and 



finally even in the negative ("magnetic") range, where (261) is quantitatively 
more accurate, the newer formula ( |263| ) still provides a description that is in 
qualitatively adequate agreement with that of its predecessor. In the ("electric") 
range x > where the formula ( |263D is most accurate, the equation of state 
relating the energy U to the tension T will take the form 

U = T + m, 2 (e 2 ^ m2 - T ^ m ' 2 - 1), (265) 

while in the ( "magnetic" ) range the corresponding relation is given by 

T = U-m, 2 (l _ c - 2(c/ - m2)/m * 2 ) . (266) 

It is evident that in the absence of background forces, knowledge of the two 
parameters m and m* is fully sufficient to characterise the dynamics of such 
a string model - indeed it is only the ratio m/m* that actually matters for 
this purpose. However one will also have to evaluate the relevant scaling factor 
tt = \f^a (which, like m and m* will depend on the parameters characterising 
the underlying bosonic field theory) if one wishes to be able to allow for the 
effect of any background electromagnetic field that may be present, since for 
this purpose it is necessary to know the value of the mathematically normalised 
coupling constant q. The latter will be directly obtainable from the value of (j 
(which in typical cases may be expected to be of the order of unity) via ( |248[ ) , on 
the assumption that the relevant value of the physical charge coupling constant 
e characterising the boson condensate is just the usual electronic charge (as 
given in Planck units by e 2 w 1/137) or a simple integral multiple thereof. 

5.8 *The self-dual transonic string model 

Although the formulae || in the preceeding section provide the simplest de- 
cription that one could reasonably desire for the purpose of describing realistic 
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superconducting strings arising from the Witten mechanism, it is instructive 
to study an even simpler model that originally turned up in a rather different 
physical context of a more artificial nature, namely that of a Kaluza Klein type 
projection of a simple Goto-Nambu model in an extended background with an 
extra space dimension in the manner first suggested by Nielsen |85| . ft can eas- 
ily be shown || that this gives rise to a particularly elegant illustration of the 
above formalism in which the relevant master function has the rather unusual 
though by no means unique property of being self dual in the sense that (for 
an appropriate choice of the scaling of the stream function ip) the dependence 
of the master function A on its argument x is exactly the same as that of the 
corresponding function A on £ their form in this case being is given in terms 
of a constant mass parameter m simply by 



A = m yj m 2 — x i A = m \J m 2 — \ ■ (267) 

This self duality property means that the equation of state relationship between 
U and T will the same form in the magnetic regime x > > % as in the electric 
regime x > > X- What is quite unique about the equation of state obtained 
for particular self dual string model is that it will have the non-dispersive - 



i.e. permanently transonic - "constant product" form (205), which can be seen 



[B9| to be expressible parametrically in terms of a dimensionless self dual state 



function i? of the kind introduced in (224) by 



U = m 2 Jl + ^— = m 2 coth -& , T = m 2 J 1 - At = m 2 tanh d , (268) 

V mr V m 2 

which gives 

c E = c L = tanhtf , (269) 



so that by (|205|) and (204), the extrinsic ("wiggle") and longitudinal ("woggle" 



or "jiggle" ) bicharacteristic vectors will coincide, having the form 



i' 



± 



coshtfw^ isinhtf^ . (270) 



The revelation H of the special transonicity property of the model char- 



acterised by ( 267 ) was sufficient to invalidate the over hasty claim []86J that 



Nielsen's elegant artifice |85 effectively represents the outcome in the "pure" 
string limit (in which effects of finite vortex thickness are neglected) of t he W it- 
ten mechanism ||. Long before the derivation j|] of the explicit formula (262) it 



was already known that in order to provide a satisfactory representation a Wit- 
ten type vacuum vortex it would be necessary to use a model (97]] of generically 
dispersive type p6[ . More specifically, it was known from the work of Peter 
[ [75| , frfjl , J77| that such a model would need to be characterised typically by 
supcrsonicity, c E > c L for spacclike and moderate timelike values of the current 
provided the ratio of the carrier mass to the Higgs mass is not too large. Now 
that the model ( ^63[ ) is available, it can be verified directly that it will indeed 
exhibit supersonicity for all spacelike values and at least for small timelike val- 
ues of the current not just when the ratio m* 2 /m 2 is small compared unity but 
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so long as it does not exceed the implausibly large value 2. On the other hand 
it can be seen that such model will always become subsonic for sufficiently large 
timclike current values (near the relaxation limit in the "electric" regime). 
It is to be remarked that despite its inability to reproduce these features, the 



use of the mathematically convenient transonic model (267), as advocated by 
Nielsen and Olescn p(| , for the purpose of describing the dynamics of a realistic 
superconducting string model was nevertheless a considerable improvement, not 
just on the use for this purpose of an unmodified Goto-Nambu model, but 
also on the more commonly used description provided by the naively linearised 
model HI , [BjjL [0 , @ , |f72|, |f74| wit h eq uation of state given by constancy 



of the trace (|208| ) which, by ( |194|) and (204), is evidently characterised by per- 
manent subsonicity, c E < c L = 1. 



The intrinsic equations of motion (223) for the transonic model can be use- 



fully be recombined J40|], |94j as an equivalent pair of d iver gence relations that 
are expressible in terms of the bicharacteristic vectors ( |270| ) , in the form 



V p [(U-T)iJ>) =-£/f p , (271) 

which shows that the "left" and "right" moving "bicharacteristic currents", 
(U — T)l + P and (U — T)£_ p will each be conserved separately in the free case, 
i.e. when f p vanishes. As an alternative presentation of the tangential force 
balance equations for this permanently transonic model, a little algebra suffices 
to show that the intrinsic equation of motion ( [226| ) can be rewritten for this 
case, in a form more closely analogous to that of the corresponding extrinsic 



equation of motion (227), as the pair of equations 



(U - T)if v V%V = « + ■ (272) 

This tangentially projected part of the dynamic equations can now be recom- 
bined [^9| with its orthogonally projected analogue (5.22), so as to give the 
complete set of force balance equations for the non dispersive permanently tran- 
scharacteristic string model (4.21) as the extremely useful pair of bicharacteristic 
propagation equations 

(U - T)£/%£/ = (g^ + £/£ T v )J v . (273) 

5.9 *Integrability and application of the transonic string 
model 



It is evident that the bicharacteristic formulation (273) of the equations of mo- 
tion for the transonic string model will take a particularly simple form when 
translated into terms of the corresponding characteristic coordinates a ± on the 
worldshcet, where the latter are defined by taking the "right moving" coordi- 
nate a+ to be constant allong "left moving" characteristic curves, and taking the 
"left moving" coordinate a~ to be constant along "right moving" characteristic 
curves, with the convention that the correspondingly parametrised bicharacter- 
istic tangent vectors t£ = dx^ fda*- should be future directed. 
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In the case of a free motion (i.e. when the force term on the right of ( (273| ) 
vanishes) in a flat background, this simplification can be used [ |B9| to obtain 
the complete solution of the dynamical equations in a very simple explicit form. 
The way this works is that with respect to Minkowski coordinates in such a 



background (273) will reduces to the simple form 

d£^/da T =0, (274) 

whose general solution is given in terms of a pair of generating curves x + {c} as 
a sum of single variable functions by the ansatz 

x» ■ = x£{a+} + x*{a-} , (275) 

which gives £^ = if? using the dot here to denote the ordinary derivatives of 
the single variable functions with respect to the corresponding characteristic 
variables. The solution ( |275[ ) generalises a result that is well known for more fa- 
miliar but degenerate Goto Nambu case, in which the tangents to the generating 
curves are required to be null, £±£ ±l i = (so that with the usual normalisa- 
tion their space projections lie on what is known as the Kibble Turok sphere 
[p8[). The only restriction in the non - degenerate case is that they should 
be non-spacelike and future directed (the corresponding projections thus lying 
anywhere in the interior, not just on the surface, of a Kibble Turok sphere) the 
unit normalisation condition fl225| ) being imposable as an option, not an obliga- 
tion, by choosing the parameter a to measure proper time allong each separate 
generating curve. 

This special property of being explicitly integrable by an ansatz of the same 
form ( |275| ) as has long been familiar for the Goto Nambu case can immediately 
be used to provide a new direct demonstration of the validity of the claim 
I originally made several years ago [^9|] to the effect that the non-dispersive 



constant product model (205) is not just a mathematical curiosity, but that it 
is usefully applicable for the purpose of providing a realistic description of the 
average motion of a "wiggly" Goto Nambu string, effectively extending to the 
tension T the concept of "renormalisation" that had previously been introduced 
for the energy density U by Allen and Shellard flQC| . 

The new justification | f)2[ summarised here is needed because my original 
argument [|39| was merely of a qualitative heuristic nature, while Vilenkin's 
mathematical confirmation |kJ was based on indirect energetic considerations, 
and was legitimately called into question f9l| on the grounds that it did not 
cover quite the most general class of "wiggles" that can be envisaged. Such an 
objection is effectively bypassed by the alternative approach presented here, but 
there would still be a paradox if allowance for more general "wiggles" really did 
produce what was alleged: thus although it was based on manifestly muddled 
reasonning and misinterpretation of preceding work, the purported demonstra- 



tion [|9l[ of higher order "deviations" from the constant product form (205) for 
the effective equation of state gave rise to a controversy that was not defini- 
tively resolved until a genuinely valid extension of Vilenkin's method ]90| was 
finally produced by Martin |93| , whose conclusion was that the constant product 
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equation of state ( |205| ) will after all provide an account that remains accurate 
for arbitrary "wiggle" perturbations of a Goto Nambu string (subject only to 
the restriction that their amplitudes should not be so large as to bring about a 
significant rate of self intersection) . 

As an alternative to the energy analysis originated by Vilcnkin and com- 



pleted by Martin, the more direct justification 92 for the use of the elastic 



string model characterised by (267) as a model for the large scale averaged be- 
haviour of a Goto-Nambu string is as follows. The argument is simply based 
on the observation that such a model implicitly underlies the diamond lattice 



discretisation that, since its original introduction by Smith and Vilenkin [ 1 1 1 , 
has been commonly employed by numerical simulators |l02| as a very convenient 
approximation scheme - of in principle unlimited accuracy - for the represen- 
tation of a Goto-Nambu string worldsheet in a flat background. As a way of 
replacing the exact continouous description by a discrete representation such as 
is necessary for numerical computation, the idea of the Smith Vilenkin method 
is simply to work with a pair of discrete sets of sampling points x£ r = x^{o~ r } 
determined by a corresponding discrete set of parameter values <j r on the gen- 
erating curves of the exact representation ( |275| ). This provides a "diamond 
lattice" of sample points given (for integral values of r and s) by 

<s = <r + ^ s > (276) 

that will automatically lie exactly on the "wiggly" Goto Nambu worldsheet 
( |275D , which is thus represented to any desired accuracy by choosing a suf- 
ficiently dense set of sampling parameter values er r on the separate "wiggly" 
null generating curves x^{a}. It is evident that the chosen set of sample points 
x£ T = x^_ {er r } on the separate "wiggly" null generators can also be considered to 
be sample points on a pair of smoothed out, and thus no longer null but timelike, 
interpolating curves that, ac cord ing to the result [|S9| demonstrated above, can 
be interpreted according to (|275| ) as generating a corresponding solution of th e 
equations of motion for an elastic string model of the kind governed by ( |205|) . 
The not so "wiggly" elastic string worldsheet constructed by this smoothing 
operation will obviously be an even better approximation to the exact "wiggly" 
Goto Nambu worldhcct than the original Smith Vilenkin lattice representation, 
which itself could already be made as accurate as desired by choosing a suffi- 
ciently high sampling resolution. No matter how far it is extrapolated to the 
future, the smoothed elastic string worldsheet generated according to ( [274|) can 
never deviate significantly from the underlying "wiggly" Goto-Nambu world- 
sheet it is designed to represent because the exact worldsheet and the smoothed 
interpolation will always coincide precisely at each point of their shared Smith 
Vilenkin lattice fl276| ). This highly satisfactory feature of providing a poten- 
tially unlimited accuracy could not be improved but would only be spoiled by 



any "deviation" from the originally proposed |89| form (205) for the effective 
equation of state. 

After thus conclusively establishing that the permanently transonic clastic 
string model characterised by the simple constant product equation of state 
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(EOq) (without any higher order corrections) provides an optimum description 
of the effect of microscopic wiggles in an underlying Goto Nambu model so long 
as self intersections remain unimportant (as was assumed in all the discussions 
[ |39"| , [ p0| , pH , Q, Q cited above), it remains to be emphasised that the 
neglect of such intersections will not be justified when the effective temperature 
p9) , (94j] of the wiggles is too high (as will presumably be the case during 
a transient period immediately following the string - forming phase transition) . 
The result of such intersections will be the formation of microscopic loops, of 
which some will subsequently be reconnected, but of which a certain fraction 
will escape. Estimation of the dissipative cooling force density that would be 
needed to allow for such losses remains a problem for future work. 



6 Symmetric configurations including rings and 
their vorton equilibrium states 

6.1 Energy-momentum flux conservation laws 

Whenever the background space time metric is invariant under the action of a 
(stationarity, axisymmetry, or other) continuous invariance group generated by 



a solution Killing of the equation (p_18j), i.e. in the notation of (106) 



££9^ = , (277) 

then any string that is isolated (i.e. not part of the boundary of an attached 
membrane) will have a corresponding momentum current (interpretable, de- 
pending on the kind of symmetry involved, as representing a flux of energy, 
angular momentum, or whatever) given by 

p M = TW . (278) 



In accordance with (120), this will satisfy a source equation of the form 

= 1^ , (279) 

which means that the corresponding flux would be strictly conserved when the 
string were not just isolated but free, i.e. if the background force were zero. 
When the string is subject to a background force of the Lorentz-Joukowsky form 



(238) that arises from background electromagnetic and Kalb-Ramond fields, 



then provided these background fields are also invariant under the symmetry 



group action generated by i.e. in the notation of (106) 



e£A ll = , ££B^ = , (280) 

it can be seen to follow that, although the physically well defined surface current 
P M will no longer be conserved by itself, it still forms part of a gauge dependent 
generalisation, say, that is strictly conserved, 

%V» = 0, (281) 
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and that is given, in terms of the gauge dependent generalisation 



T"„ = T"„ + + B^S„ P (282) 
of the surface stress momentum energy density tensor, by 

= = P p + (qcfAv + hS^Bup)^ . (283) 

6.2 Bernoulli constants for a symmetric string configura- 
tion 

Let us now restrict attention to cases in which string configuration itself shares 
the background symmetry under consideration. This category includes the as- 
trophysically interesting case of circular configurations in an axisymmetric back- 
ground - which is of particular interest in the context of vorton formation in a 
flat background - and of course it also includes many familiar terrestrial exam- 
ples of static string configurations in a stationary background - one of the most 
obvious being that of the overhead telephone and power cables that festoon our 
environment. Any string configuration that is symmetric in this sense will be 
characterised by the condition 

_L^r = , (284) 

meaning that the relevant symmetry generator is tangential to the worldsheet, 
and the corresponding Lie invariance condition on its surface stress momentum 
energy density tensor will have the form 

l p XljF vp = 2T^V f ,£ p) . (285) 

Under such conditions, as well as the ordinary momentum flux P' 1 , what may 
be termed the adjoint momentum flux, 

fp> = ^tjy* j tjy* = E vp E^T p a , (286) 



will also obey an equation of the form (279), i.e. 



%P" = , (287) 

and hence would also be conserved if the string were free. In the presence of 
an electromagnetic or Kalb Ramond background, it can be seen that, like the 
ordinary momentum flux, this adjoint momentum flux has a gauge dependent 
extension, 

fpM = t \Y u p , *T V » = £ vp £ ii(J r p a , (288) 

that will share with the generalised momentum flux V p (that would be conserved 
even if the string did not share the symmetry of the background) the property 
of obeying a strict surface current conservation law, namely 

% f P'' = . (289) 
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The group invariance conditions 



r v„ P/1 + P „ v = o , r v„p M + p„v M r = o , (290) 



that are the analogues of (285) for the separate mutually dual pair of internal 
momentum covectors p^ and associated with the internal current within 
the string, can be rewritten, with the aid of the corresponding electromagnetic 



background invariance condition fl28Cf) , in the form 

£ p £^ (V v p„ + |f„„) = £ pv W v B , iPE^VvPr = £ pu V v p , (291) 

in terms of a pair of scalars B and /3, that can be considered as generalisations of 
the well known Bernoulli constant in stationary classical perfect fluid dynamics, 
and that are defined by 

B = P + qA^ , = p^ , 13 = p^ . (292) 

If V 1 is timelike, the corresponding symmetry will be interprctable as station- 



arity, while the more restrictive case [102] in which the string is actually static 



(in the sense that there is no transverse current component relative to the back- 
ground rest frame determined by will be given by the condition that the 
second Bernoulli constant, (3 should vanish. 



Comparing the conditions ([291]) with (244) and (246) it can be seen that 
the internal dynamical equations are equivalent in this group invariant case to 
the corresponding pair of Bernoulli type conservation laws to the effect that B 
and (3 (but, unless the electromagnetic field is absent, not (3) should both be 
constant over the worldsheet. This observation allows the problem of solving 
the dynamical for a symmetric configuration of the kind of ( "perfectly elastic" 
, i.e. barotropic) string model under consideration to one of solving just the 
extrinsic equations governing the location of the worldsheet. 

6.3 Generating tangent vector field for symmetric solu- 
tions 



A recent investigation [104] based on the systematic use of variational meth- 
ods in the restricted case for which the Kalb-Ramond coupling was absent has 
drawn attention to the interest of extrapolating the Bernoulli constants out- 
side the supporting worldsheet as a pair of scalar fields defined over the entire 
background spacetime by the uniformity conditions 

V^B = , V M /3 = , (293) 

and formulating the problem in terms a certain preferred generating vector 
field that is constructed in such a manner as to be everywhere tangent to the 
worldsheet according to the specification 

x » = S^^p^ = te vp T pp . (294) 



64 



Using (251) and (292), it can be seen that the generating vector defined in 
this way will be expressible in terms of the variables introduced in Section |^ as 

1" = ^ + ^, (295) 

X X 

while the Killing vector itself will be expressible in analogous form by 

F = — d 1 + -J— . (296) 
A-A A-A 

By further use of the formulae in Section it can be seen that the squared 
amplitude of the latter (which in the case of stationary symmetry can be inter- 
preted as an effective gravitational potential field) will be given by an expression 
of the form 

B 2 B 2 

e = m u = -— - . (297) 

XX 

To obtain a solution of the dynamical equations governing a symmetric string 
configuration of he kind under consideration it will evidently be sufficient to 
obtain an integral trajectory of the generating tangent vector X^, since when 
this is available the complete worldshhct can then be swept out in a trivial 
manner by dragging this trajectory along under the symmetry action generated 
by . A very useful approach to the problem of finding sych a tange nt tra jectory 



is based on a procedure that exploits of the prccceding equation ( 297 ) in the 
following manner. 



In view of (292), and of th e st ate functional relationship (|24fj) between x 
and its dual, the equation (|297j ) can be solved for any particular choice of 
the constant "tuning parameters" B and (3 to determine the internal variables 
X and x as functions of the (geometrically determined) potential £ 2 and the 
(electric type) potential t v A v , and hence by implication as scalar fields over the 
entire background space, not just on the worldsheet where they were originally 
defined. 



Again using the formulae in Section |L it can be seen from (295) and ( [296| 



that the contraction of the worldsheet generating vector A" M with the Killing 
vector field will be given simply by 

PXp = 0/3 , (298) 

and that the squared amplitude of will given by an expression of the form 

X 2 = X»X„ = + . (299) 

X X 

By the latter equation, X^X^ is also implicitly defined as a function of t 2 
and i v 'A v , so that it too can be considered defined scalar field, not just on the 
worldsheet but also over the background as a whole. The gradient of the field 
constructed in this way will be given by the expression 

V^X 2 = -AA V„£ 2 + 2qj3(A - ~h)^V»{l v A v ) . (300) 

X 
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Previous experience [104] with the case in which only the electromagnetic but 
not the Magnus force contribution is present suggests the interest of formulating 
the problem in terms of the propagation of the special generating vector AT M , 



which (using the formulae ( |210| ) and (251) of the previous section) can be seen 
from (|295|) and (E9q) to be given by 



X^X^ = A A t u V u l p + X v Z vp l p ] p , 



(301) 



where is the background force as given by (237). It can now be seen that 
the two preceeding equations can be combined to give the equation of motion 
for the worldshcct generating vector X^ in the very elegant and convenient final 
form 

X V V V X„ - \v p X 2 = T plJ X v , (302) 



in terms of a pseudo Maxwellian field given by 



(303) 



where A p is a gauge dependent pseudo-Maxwellian potential covector given by 



An = qPA^ + hB 



(304) 



If is timelike so that the corres pond ing symmetry is interpretable as sta- 
tionarity, then the equation of motion (302) will be interpretable as the condition 
for the string to be in equilibrium with the given values of the constant "tuning" 
parameters B and (3, (of which, as remarked above, the latter, ( 3, wi ll vanish in 
the case of an equilibrium that is not just stationary but static, [ 103 ) . The new 
result here is that the Jou kowski type "lift" force (which was not allowed for in 
the previous analysis [ 104 1 ) due to the Magnus effect on the string, as it "flies" 
(like an aerofoil) through the background medium represented by the current 
3-form N^p, has just the same form as an extra Lorentz type electromagnetic 
(indeed in the stationary case purely magnetic) force contribution. 



6.4 Hamiltonian for world sheet generator 

When the genuine electromagnetic background coupling, and the similarly act- 
ing Kalb Ramond coupling are both absent, then (as pointed out previously 
[104 1) the equation of motion ( |302|) (that is the equilibrium condition for "steady 
flight" in the stationary case) is just a simple geodesic equation with respect to, 
not the actual background spacetime metric g pv , but the conformally modified 
metric X 2 g liv , with the conformal factor X 2 determined as a field over the back- 



ground by ( |299| ) in conjunction with (293) and (297). Eve n wh en the Lorentz 
and Joukowski force contributions are present, the equation (302) governing the 
propagation of the world sheet generator X^ retains a particularly convenient 
Hamiltonian form, given by 



A"' 



dx^ dH 



dT 



dn L 



dn p 

dr 



dH 



(305) 
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for the quadratic Hamiltonian function 

h = \ 9 ^{n lx -A lx ){n v -A v ) 



:X 2 



(306) 



subject to a restraint fixing the (gcncrically non-afhnc) paramctrisation r of the 
trajectory by the condition that the numerical value of the Hamiltonian (which 
will automatically be a constant of the motion) should vanish. 



H = 



(307) 



together with a further momentum restraint, determining the (automati cally 
conserved) relative transport rate II ^ in accordance with the relation ( [298|) 
by the condition 

Upl* = B(3 . (308) 
The Hamiltonian momentum covcctor itself can be evaluated as 



A — tf T p 



(309) 



In terms of the original conserved generalised momentum fluxes given by ( [283| ) 
and (288, and the gradients of the scalar (stream function and phase) potentials 
introduced in (230) and (247), this Hamiltonian covector will be given by 



n, 



A ¥ 



f 



\t> v 



B\,,,r 



(310) 



(It is to be noted that there were omissions due to transcription errors in the 
analogues of the latter formulae in the preceeding version (TJ of these notes.) 

The advantage of a Hamiltonian formulation is that it allows the problem 
to be dealt with by obtaining the momentum covector in the form II ^ = V^S 
from a solution of the corresponding Hamilton Jacobi equation, which in this 
case will take the form 



gT ( V M 5 - A») { V„S - A v ) = X 2 



(311) 



with X given by (299) via (297), while the restraint (308) gives the condition 



rv„s = QB 



(312) 



It is to be remarked that the square root X of X 2 may be purely imaginary 
since the generating vector may be timelike. Instead of working with X it has 
been found convenient |l06fl to work with a related field, T, that will be real if the 
Killing vector is spacelike so that I is real, and that is given explicitly by either 
of the mutually dual (by (23G) and (297) identically equivalent) alternative 
formulae 

p 2 ic ~ p 2 ic 



r 



A 



M. 



(313) 



Considered as a function of the field t defined by (297) the derivative of T will 
be given by 

(314) 



d( 



(2 
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which can be seen to be equivalent to 



° 2 ll = 



It can be seen that the field X 2 in (311) will be expressible in terms of this new 
self-dual field Y = T in the form 

B 2 B 2 

X 2 = ^--Y 2 . (316) 
6.5 Free evolution of circular loops 

A particularly simple and important application of the formalism that has just 



been presented is to case [ 105 1 , [106 of a small circular string loop undergoing 
motion that is free in the strong sense, meaning that there are no external forces 
of axionic, electromagnetic or even gravitational origin, i.e for which B M „ and 
can be taken to vanish, and for which can be taken to be flat, which 
will be a realistic approximation under a wide rang of circumstances, including 
cosmological applications in which the loop scale is small compared with the 
relevant Hubble length scale. There will be no loss of generality in taking such 
a circular loop to lie in the equatorial plane of a spherical coordinate system 
with origin at its center of symmetry, in terms of which the metric will have the 
standard form 

g^dx^dx" = dr 2 + r 2 (d6 2 + sin 2 0d(f> 2 ) - dt 2 . (317) 

This means that the loop lies in the hyperplane 6 = 7r/2, with a worldsheet 
whose internal coordinates could be taken to be t and <f>. The extrinsic location 
of such a worldsheet will be specifiable just by giving its radius r as a function 
t. The Killing vector £ M generating the symmetry of the world sheet can be 
conveniently taken to be given by 

r ^= 2 4' (318) 



so that its magnitude I as given by fl297| ) will be directly identifiable with the 
circumference 

i = 2nr (319) 

of the loop. The flux of angular momentum as defined with respect to the 
usual angular momentum generator d/d(f> will thus be given by the proportion- 
ality relation 

P M = 2ttJ^ , (320) 



where P M is the momentum curre nt (278 ) associated with £^ (which in this 



case is the same as that given by (283), since we are neglecting the effects of 



electromagnetic and axionic fields). The stationarity of the background (though 
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not of the worldsheet) with respect to the Killing vector that generates time 
translations according to the specification 



k" 



d 



d_ 

dt' 



(321) 



ensures that as well as the conserved angular momentum flux ( p20| ) there will 
also be a non-trivially conserved mass-energy flux vector that is specified by the 
corresponding analogue of (283) as 



(322) 



For each of these independent conserved surface fluxes there will be a corre- 
sponding conserved constant characterising the loop, namely its total angular 
momentum 

J = j dx» E^J V , (323) 



and its total mass 



M = <j> dx* 1 £ ^M" 



(324) 



As well as these two quantities, whose conservation depends on the pos- 
tulated summctry, there will be two other integrals that will be conserved in 
any case, namely th e wi nding num ber N and the particle quantum number Z 
that are defined by ( 256 ) and (257), whi ch ca n be seen to be determined in the 
present case by the Bernoulli constants ( 292 ) , in terms of which they are given 

by 

Z = 2ir$N = B . (325) 

The angular momentum and mass integrals can also be expressible by purely 
local formulae, in terms of the the momentum covector H M for the worldsheet 
generator. Since the effective gauge field vanishes in the present application 
this momentum covector will be given directly by the formula (p94T) as 



n n = ts vp T" 



f ■ 



(326) 



from which it can be seen that the corresponding total angular momentum and 
mass constants will be given by 



PJ7„=2ttJ, 



-M . 



(327) 



It is to be noticed that only three of the f our cons erve d integrals J, M, Z, N 
are actually independent: by substituting ( |236[) in ( 326 ) it can be seen that the 
first one is just the product of the last two, i.e. 



J = ZN . 



(328) 



(It it to be remarked that yet another way of obtaining a conserved integral 
would to take the adjoint analogue of the angular momentum that is con- 
structed from t as given by (J28q) instead of from as given by (p7q) 
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according to the specification 2n^J = l^E^v ^ P v , but the result in this case is 
not independent: it can be seen that the angular momentum is self adjoint in 
the sense that t J = J.) 



The mass-energy constant (327) is intcrpretable in the present case as the 
rate of variation of the ordinary coordinate time t with respect to the Hamilto- 
nian parameter r introduced in (|305[) : 



(329) 



Using a dot to denote differentiation with respect to the coordinate time t, the 
radial momentum component II \ will correspondingly be given by 



dr 



Mr. 



In terms of this quantity the Hamiltonian (306) will be expressible as 



M - X z 



(330) 



(331) 



I t is at thi s po int that the identity (316) is useful, since it can be seen from 
( |325|) and (|328| ) that it provides the convenient simplification 



J 2 2 



(332) 



The normalisation constraint (307) to the effect that the Hamiltonian should 
vanish thus provides the required equation of motion for the loop radius r in 
the conveniently first integrated form 



2S2 



M - T 1 



(333) 



where T is the function specified - independently of M - by the Bernoulli 
constants (i.e. by Z and TV) according to the formula (313), whose explicit form 
depends of course on the equation of state. 

For the case of a string model representing a Witten type superconducting 
vortex an app ropriate approximation for the equation of state is provided by the 
formula ( |263| ) whose application to this problem has recently been the subject 
of detailed investigation ||l06f| (following a more specialised study of the zero an- 
gular momentum case for which Z or N vanishes so that collapse is unavoidable 
[107 1): the function T has no maximum, but for moderate values of the ratio of 
Z to N it has a finite minimum, and in that case, depending on the particular 
values of the three independent constants Z, N, and M the loop may either 
oscillate smoothly between a finite maximum and minimum value of the radius 
or else it may reach a critical minimum radius at which the model breaks down 
due to overconcentration of the current which results in a local instability: in 
the case of a timelike current such an instability arises because the "wiggle" 
speed c E given by (194) tends to zero, whereas in the case of a spacelike current 
it arises because the longitudinal characteristic speed c L given by (204) tends 
to zero. 
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6.6 Vorton equilibrium states 

It can be seen from the final equation fl333p of the Subsection |6.5| that the circular 
loop will have a stationary equilibrium state characterised by the minimum 
admissible value of M for given values of Z and TV at the radius r for which 
T takes its minimum value. This happens for M = T where dT/dr vanishes, 



which according to (297) and ( pl4[ ) occurs where 

(3 2 Ax + P 2 A X = . (334) 



According to (315) this is equivalent to 

T^i^ = (335) 

which is interpretable as the condition that the stress in the direction of the 
axisymmctry generator should vanish. 

Such ccntrifugally supported ring states are simple examples of the phe- 
nomenon for which Davis and Shellard |8lJ , |82| have introduced the term "vor- 
ton" . The possibility of forming "vortons" , meaning ccntrifugally supported but 
not necessarily circular equilibrium states of cosmic string loops is of consider- 
able potential importance for cosmology, but it was largely overlooked during 
the early years of cosmic string theory. 

The formalism embodied in the generator equation ( |302| ) and the corre- 
sponding Hamiltonian ( |306| ) can be applied to the general study of stationary 
equilibrium states by replacing the axisymmetry vector £ to which it was applied 
in the preceedmg subsection by the Killing vector fc M that generates the station- 



ary symmetry of the background, which, when substituted for in (294), will 
determine a corresponding world sheet tangent vector 

a" = k"£ up T p ^ , (336) 

instead of the vector X 1 *. 

Such an approach provides results of considerable mathematical (as opposed 
to physical) interest when applied to the case of stationary strings in a black 



hole background. Generalising results obtai ned previously 108 , |109| for the 



Goto Nambu limit case, it has been shown [ 104 1 that for a string model of the 
non-dispersive permanently trasonic type with the constant product equation 
of state ( [205| ) that is governed by the Lagrangian (p67|) (as obtained § both 



from the Nielsen dimensional reduction mechanism and also [|39| , as explained in 
Section 5, from the more physically realistic "wiggly" string approximation) the 
stationary Hamilton Jacobi equation is exactly soluble by separation of variables 
in a Kerr black hole spacetimc, not just of the ordinary asymptotically flat kind 



but even of the generalised asymptotically De Sitter kind [110], 1 1 1 1 ] . Except in 
the Schwarschild-De Sitter limit, where it could of course have been predicted 
as a consequence of spherical symmetry, this separability property still seems 
rather miraculous, reflecting a "hidden symmetry" of the Kerr background that 
is still by no means well understood. The newly discovered separability property 



[ 1 04 1 is not just an automatic consequence of the simpler, though when first 
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discovered already surprising, property of separability for the ordinary geodesic 



equation [112 but depends on a more restrictive requirement of the kind needed 



for the more delicate separability property of the scalar wave equation [ 1 1 3 1 



(It is however more robust than the separability properties that have turned 



out to hold for higher spin bosonic |114| and fcrmionic [115], ]116[, [117] wave 



equations, and other related systems JL 18| , [119].) 

From a physical point of view the mathematically interesting configurations 
discussed in the preceeding paragraph are rather artificial. What is much more 
important from a cosmological point of view is the the equilibrium of small 
closed string loops in the mathematically relatively trivial case for which the 
background gravitational, electromagnetic, and axionic forces are negligible. In 
the absence of electromagnetic, or axionic forces, the relevant generator equa- 



tion, which is obtained (by substituting for in (302) with the source term 
on the right set to zero) will have the form 

ff"V„<7 M - \^{a v a v ) = • (337) 

It is evident that the worldsheet generators obtained as solutions of this equa- 
tion in a Minkowski background will all just be straight lines for the case of 
a stationarity Killing vector k, since in this case (unlike that of axisymmctry) 
k v k v and hence also a v a v will be constant. Since in this case the Killing vector 
trajectories are also straight, it might at first seem to follow that the worldsheet 
of a stationary string in an empty Minkowski background would necessarily be 
flat. This conclusion would exclude the possibility, when on background field 
is present, of any closed loop equilibrium states - such as the circular con- 



figurations characterised by (334). In fact cosmologists did indeed (albeit for 
other reasons) entirely overlook the possibility that such states might exist until 
the comparitively recent publication of an epoch making paper by Davis and 
Shellard flsi| provided the first counterexamples (the only previously considered 
equilibrium states |12(| , 170] , Q , Wa having been based on a magnetic support 



mechanism that was was finally judged to be too feeble to be effective except 



[122 1 as a minor correction). 

The loophole in the deduction that if the trajectories generated by cr M and by 
fc M are both straight then the worldsheet must be flat is that it is implicitly based 
on the assumption that the two kinds of trajectories cross each other transversly. 
However there will be no restriction on the curvature in the transverse direction 
in the critical case for which the two kinds of trajectory coincide, i.e. for which 
a* 1 and fc M are parallel. The condition for criticality in this sense is expressible 
as 

a^E ^k" = , (338) 

which can be seen from the definition (|336| ) to be interpretable as meaning that 
the extrinsic characteristic equation ( |133| ) is satisfied by the tangent covector 
Xfi = Eitvk", which is proportional to l^. This is equivalent to the condition 
that the Killing vector fc M itself should be bicharacteristic, in the sense of being 
directed allong the propagation direction of extrinsic perturbations of the world 
sheet. The criticality condition (|338|) - which can be seen to be equivalent to 
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( |334D - is thus interpretable as a condition of characteristic flow. It means that 
the "running velocity" , v say, of relative motion of the intrinsicly preferred rest 
frame of the string (as determined by whichever of c M and is timelike) relative 
to the background frame specified by the (in this case necessarily timelike) 
Killing vector fc M is the same as the extrinsic propagation velocity c E given by 



®, i-e. 

T 



,2 



V 



(339) 

In the presence of generic gravitational, electromagnetic and Kalb Ramond 
forces, the criticality condition ( |33§| ) can be satisfied at particular positions, such 
as where there is a transition from a subcharacteristic running velocity, v < c E 
to a supercharacteristic running velocity v > c E (as will occur for instance on a 
string in a steady state of radial flow into a Schwarzschild black hole) or where 
there is a cusp, with v = c E but with subcharacteristic flow v < c E on both 
sides. However in view of the constancy of the Bernoulli "tuning" parameters 
B and (3, the absence of any back ground field will always allow, and generically 
(the special integrable case (p68|) being an exception) will ensure, uniformity of 



the state of the string, so that the transcharacteristic flow condition (338) can 
be satisfied throughout its length. The space configuration of such a uniformly 
transcharacteristic steady string state can have arbitrarily variable curvature, 
and so is compatible with a closed loop topology. 

The value of the mass M at which the equilibrium condition ( |339| ) is attained 
will of course be proportional to the circumference £ of the vorton: it can easily 
be seen that its mass to length ratio will be given by 

M 

— = U + T = -A- A. (340) 

For given values of the conserved quantum numbers N and Z, the absolute 
scale of the vorton will be given by the formula for its length which takes the 
comparatively simple form]83|] 

, 2= ^ra = «. (341) 



In the "chiral" case when the current is approximately null, and more g ener ally 



in the magnetic case for which it is spacclike, one expects (see Subsection 5.7 ) the 
energy density and tension to be given by U w T r* m 2 , where m is the relevant 
Higgs mass scale. On this basis the two preceeding formulae are sufficient by 
themselves to provide reasonably accurate estimates of the vorton length and 
circumference. However if the current is timelike, U may become relatively large 
and T very small compared with m 2 . To estimate their values in a case of this 
more general kind, it is necessary to solve the (non linear) equilibrium equation 
( |334|) which (since our assumption that electromagnetic effects are unimportant 
implies (3 = B) determines the local state as a function of the ratio of the particle 
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number Z to the winding number N by a relation of the form 



(342) 



6.7 * Allowance for the electromagnetic "spring" effect. 

Instead of working with the generator tangent vector <r M that is given by the 
(not so obvious) specification ( |336D , another way of deriving the conclusions 
of Subsection 6-6, and in particular the centrifugal equilibrium condition ( |339D , 
is to use an approach whereby the quantities in the fundamental dynamical 
equation ( |122| ) are directly evaluated in terms of the worldsheet tangent vector 
e p that is given modulo a choice of sign by the (more obvious) specification that 
it be orthogonal to the Killing vector fc M with unit normalisation, 

e v k v = , e v e v = 1 . (343) 

Since the Killing vector generator of ordinary stationary smmetry in flat space 
is covariantly constant and can be taken to have unit normalisation 

V M Av = 0, k u k v = -l, (344) 

the first fundamental of the stationary string worldsheet will be expressible by 

if = -k p k v + e"e" , (345) 

so that by (|lj) the second fundamental tensor will be obtainable simply as 

K^ f = e„e v K" , (346) 

where 

K p = e a V a e p . (347) 
Under these conditions, the preferred orthonormal diad introduced in Section 



5.2 will be expressible in the form 

u" = (1 -v 2 y 1/2 (k» + ve p ) v^ = (l-v 2 )- 1 / 2 (e tl + vk^), (348) 

where v is the longitudinal "running velocity" characterising the rate motion of 
relative to the stationary background of the internally preferred rest frame of 
the string that is determined generically by whichever of c M and c M is timelike. 
Using the expression ( |99"| ) for the surface stress-energy tensor, the left hand side 
of the fundamental dynamical equation (|122]) can thus be evaluated as 



— fUv 2 -T\ 

In the absence of external - electromagnetic or Magnus type - forces, the 
condition for equilibrium is that the vector given by (|345|) should vanish. This 
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requirement can be satisfied in one or other of two alternative manners. One 
way is for the configuration is to be straight, so that K p = 0, in which case 
any value (including the static case v = 0) is possible for v. However a straight 
configuration is possible only if the string is infinite or else terminates at fixed 
endpoints. For a closed string loop in flat space the curvature vector K p must 
necessarily be non-zero at least along part of the configuration, and in this case 
it is evident that the only way for the equilbrium condition to be satisfied is for 



the running velocity to satisfy the transcharacteristic condition ( 339 ), in which 



case there will be no restriction at all on the curvature vector K p . Since by 



(194), the centrifugal equilibrium condition (33E) is interpretable as meaning 
that the running velocity must be equal to the velocity of transverse ( "wiggle" ) 
type perturbations, it can be seen to follow that relatively backward moving 
perturbations will have their flow speed exactly canceled out, i.e. they will just 
be static deformations, which explains why, in this transcharacteristic case the 
curvature K p is unrestricted. 



The formula (349) will of course still be applicable in more general cases, 



for which the dynamical equation (122) has a non vanishing right hand side ex- 



pressing the effect of an electromagnetic or Magnus type force (the latter being 
responsible for the equilibrium of an ordinary "smoke ring" , which is the most 
familiar laboratory example of a vorton type equilibrium state). In a cosmo- 
logical context particular interest attaches to the case of electromagnetic self 
interaction arising in the case of a current with a non negligible electromag- 
netic charge coupling. Unlike the effect of a Magnus force or a purely external 
magnetic field which are straightforwardly tractable within the framework de- 
veloped here, allowance for electromagnetic self interaction raises an awkward 
and delicate problem because the self induced field is singular on the worldsheet, 
so that its treatment (as in the more familiar case of a particle self interaction) 
requires the use of some sort of renormalisation procedure, typically involving a 
cut off length representing the finite thickness of the worldsheet in a microscopic 
description. 

No satisfactory treatment of electromagnetic (or gravitational) self inetaction 
has yet been developed for the general, dynamically evolving case, but for the 
very special case of a stationary circular string loop an adequate analysis of the 
corrections arising from "spring" effect due to electromagnetic self interaction 
has been provided by Peter ||34|| . This treatment involves the introduction of a 
cut off length, r* , say that represents the effective thickness of the microscopic 
current destribution in the string and that is expected to be of the order of 
magnitude of the Compton wavelength associated with the Witten carrier mass 



scale m* that was introduced in Section 5.7, which (in Planck units) means 
r» ~ 1/m*. In the limit when this radius goes to zero the electromagnetic field 
is logarithmically divergent, its effective mean value within the string worldsheet 
being given in terms of the surface current vector j p by an expression of the 
form 

= -21n{-^}(/i [/i fc H fc a + K {ll e v] e a )f , (350) 
where r is the the radius of the circular configuration under consideration, which 
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will be given in terms of the corresponding curvature vector by 

K p K n = -4 . 



(351) 



In terms of the number den sity cu rren t vector d 1 and the associated momentum 
vector p M as given by (230) and (234) the corresponding electromagnetic force 
vector will be given by 



K h 



(352) 



where q 2 is the effective charge coupling constant which will be given in terms 
of the charge per fundamental particle, e 2 ~ 1 /137 by q 2 = § 2 e 2 where jj is the 
normalisation factor introduced in Subsection 5.6, which is expected to be of 
order unity. 

It is convenient to rewrite this in terms of the modulus \\\ of the squared 
current magnitude x = c^c^ of the current vector, which, according to the 
results described in Subsection 5.7, is expected to be subject to a saturation 



limit given in terms of the Witten mass scale to* by x 
resulting expression 



< 



TO* 



F i l 



K, 



Setting the 



(353) 



equa l to the term (349) in accordance with the fundamental dynamical equation 
( 122 ), it can be seen that the curvature covector K p cancels out again, as in 
the electromagnetically uncoupled case. The non-trivial equilibrium c ondi tion 
that remains is Peter's generalisation |3 of the preceeding formula (339) for 
the running velocity v, which is expressible in the form 



T 



U 



(354) 



where the correction term allowing for the electromagnetic self interaction is 
given by 

(355) 
7\ 



In any realistic cosmic string model of the kind described in Subsection 
it is to be anticipated that we shall always have 



X 



and one also expects to have 



< 



* 2 <TO 2 <[/., 



1/137. 



(356) 



(357) 



More specifically, in the "magnetic" range, \ > 0, for any such model, it is to 
be expected that the tension T will remain comparable with the energy density 
U , which implies that the condition 



•C T : 



U 



(358) 
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will be satisfied - unless the logarithmic factor is enormous, which is only con- 
ceivable for implausibly large values of the radius (too large for the ring to be 
treated as an isolated system unaffected by external matter). The electromag- 
netic force contribution will therefore be relatively negligible, not only in the 
"chiral" case, i.e. when current is approximately null, but also in the "mag- 
netic" case considered in early work, i.e. when the current is spacclike. This 
means that the original notion jn| , |7lJ , [|72| of non rotating magnetostatically 
supported "cosmic spring" configurations was unrealistic ||77| . 

On the other hand however, as Peter has pointed out p4), if the definition 
of a "cosmic spring" is extended to include non-rotating configurations with 
electrostatic as opposed to magnetostatic support, then the preceeding objection 
no longer applies. This is because, unlike the spacelike case, the tension can 
become relatively small, T <C U, in the case of a timelike current for a cosmic 



sting model of the kind described in Subsection p.7\ While the existence of 
exactly non rotating equilibrium states which, according to (354), would be 
characterised by the very special electrostatic support condition 

T = e 2 . (359) 

will in principle be possible, nevertheless the formation of such "spring" type 
vortons would presumably be statistically rare under natural conditions. There 
will however be a more reasonable chance for vortons to be formed in the ex- 
tended range of nearby slowly rotating states in which, if not actually dominant, 
electrostatic support does provide a contribution that may be non negligible 
compared with the purely centrifugal support mechanism considered Subsection 



6.6. Although one would expect that such very slowly rotating states would ini- 
tially represent only a small fraction of the total vorton population, they might 
come to be of dominant importance in the long run, due to the greater vulncra- 
blity of supersonically rotating states to secular and even dynamical instability 

S d, & 

6.8 *Vortons in cosmology 

The question of closed loop equilibrium states did not arise in the earliest stud- 



ies 1 121 of cosmic strings, which were restricted to the Goto-Nambu model 
whose bicharacteristics are always null and so can never be aligned with the 
timclikc Killing vector generating a stationary symmetry. However in a generic 
string model (3(| for which the 2-dimensional longitudinal Lorentz of the inter- 
nal structure is broken by a current, whether of the neutral kind exemplified 
in an ordinary violin string or the electromagnetic kind exemplified [p7| , |f75f , 
[ |T6| , J7?j by Witten's superconducting vacuum vortex model, the bicharacteris- 
tic directions will generically be timelike [^lj (spacelike bicharacteristics being 
forbidden by the requirement of causality) so there will be no obstacle to their 
alignment with a timclikc Killing vector in accordance with the criticality condi- 
tion (|338|), i.e. to having a running velocity given by v = c E . The simple "toy" 



complex scalar field model on which the pioneering cosmic string studies [|J wore 



77 



based had longitudinally Lorentz invariant vortex defects (of "local" or "global" 
type depending on the presence or absence of coupling to a gauge field) that 
were describable at a microscopic level by string models (with Kalb Ramond 
coupling in the "global" case) that were indeed of the special Goto-Nambu type. 
However the extra degrees of field freedom (starting with the additional scalar 
field introduced by Witten in his original superconducting example || ) that are 
neede d in su cces sively more realistic models |122[ , [|123| make it increasingly dif- 
ficult [ |124| , [125 1 - though it may still be possible |126[| - to avoid the formation 
of internal structure breaking the longitudinal Lorentz invariance and reducing 
the extrinsic characteristic velocity to the subluminal range c E < 1 at which 
stationary equilibrium with the critical running velocity v = c E < 1 becomes 
possible. 

The cosmological significance of this is that whereas Goto Nambu string 
loops cannot ultimately avoid gravitationally or otherwise radiating away all 
their energy [ 127 1 , 128], [ 129 [ , since they have no equilibrium states into which 
they might settle down, on the other hand more general kinds of strings, whose 
occurrence would now seem at least as plausible 1 122 , [ 123 1 , 1 124 , Hl2E 



leave a relic distribution of stationary loop configurations that may survive in- 
definitely @, @, pi, pi], p2|, pi, pl|. One would expect such con- 



figurations to be those that minimise the energy for given values of the relevant 
globally conserved quantities, of which there might be a considerable number in 
the more complicated multiply conducting models |)4| that might be considered, 
but of which there are only a single pair in the "barotropic" type string con- 



sidered here, namely the quantities Z and N defined in Subsection 5.6, which 
are proportional to the stream function winding number § dip and the phase 
winding number § d^>, whose respective constancy results from the conservation 
of the mutually dual pair of currents c M and c M . If the equilibrium is predomi- 
nantly due to a single current with corresponding particle and winding numbers 



Z and N, then according to the formulae (340) and (341), the corresponding 



vorton mass M and length I will have an order of magnitude given in terms of 
the relevant Higg's mass scale to by rough estimates of the form 



M w \ZN\ 1/2 m 



<k \ZN\ 



(360) 



whenever the current is spacelike (the "magnetic" case), or approximately null 
(the "crural" case), so that one has T U « to 2 . The exceptions - for which 
considerably higher values may be obtained - are very slowly rotating low ten- 
sion states, with a strong timelike current, which according to (342) will occur 
only for Z 2 3> N 2 : this would presumably be comparatively rare on the basis 
of random statistical fluctuations, though it might ultimately be favoured by 
selection if, as is not implausible, the initially more common (supersonically 
rotating) states with Z 2 comparable to TV 2 , turn out to be relatively unstable 
in the long run. 

The potential cosmological importance of such a distribution of stationary 
relic loops, referred to as "vortons" , was first pointed out by Davis and Shcllard 



1 13C ] , who emphasised that in the case of the "heavyweight" cosmic 
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strings whose formation during G.U.T. symmetry breaking had been postulated 
to account for galaxy formation, the density of the ensuing "vortons" would be 
more that sufficient to give rise to a catastrophic cosmological mass excess (of 
the kind first envisaged as arising from the formation of monopoles) even if they 
were formed with very low efficiency. This prediction was implicitly based on the 
assumption that the Witten carrier mass scale, represented by the parameter 
m* in ( |263| ), was of the same order of magnitude as the relevant Kibble mass 
scale, to, determined by the Higgs field, as would be the case if both the string 
formation and the current condensation were part of the same G.U.T. symmetry 
breaking phase transition, for which the corresponding gravitational coupling 
constant would be given by v Gm 2 w 1CP 3 . 

Subject to the foregoing assumption - i.e. that the Witten mass scale to* 
was not much smaller than the Kibble mass scale to - early quantitative es- 



timates suggested [131], [132 that, to avoid a cosmological mass excess if the 



vortons were formed with relatively high efficiency, the relevant phase transition 
energy scale could not have been much greater than a limit given very roughly 
by v Gm 2 ~ 10~ 13 , which is not so very much larger than the scale of elec- 
troweak unification. This limit has however been subject to upward revision by 
more recent and detailed investigations, including allowance for more realistic - 



i.e. lower - estimates of the efficiency of vorton formation ]133[, 1 134 ; the limit 



is further relaxed if one drops the supposition that the Witten mass to* should 
be close to the Kibble mass scale to, so that provided VGto* 2 > 10~ 10 string 
formation at the G.U.T. scale vGro 2 w 10~ 3 is cosmologically admissible after 
all even if the ensuing vortons are stable enough to survive to the present epoch. 
If the vortons are only stable enough to last for a few minutes, the condition that 
the standard nucleosynthesis pr ocess should not have been disrupted provides 
a somewhat weaker limit given pi roughly by y/Gm* 2 > 1CT 12 . It follows 
from this work that if the current condensation does not occur until the epoch 
of electroweak symmetry breaking then the ensuing vortons would contribute 
at most a very small fraction of the present cosmological closure density. This 
conclusion applies in particular to scenarios in which the relevant strings are 
themselves formed during electroweak symmetry breaking - which would hap- 
pen only in a non-standard (e.g. supersymmetric) theory, since the standard 
Glashow Weinberg Salam version, does not give rise to stable string like vortex 



defects [135 . However even if their density is low by cosmological standards it 



does not necessarily follow that the vortons would be undetectable: for exam- 
ple it has recently been suggested [ 136 1 that they might account for otherwise 
inexplicable cosmic ray observations. 

The potential importance of such cosmological and astrophysical effects pro- 
vides the motivation for more thorough investigation of equilibrium states that 
may be involved, a particularly important question being that of their stability. 
Prior to the derivation of the general symmetric string generator equation ( |302| ), 
the only closed loop equilibrium states to have been c onsi dered were the circular 
kind to which the use by Davis and Shellard J8l[], [130|of the term "vorton" 
was originally restricted. The first general investigation M, [§31 of such circular 
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"cosmic ring" states showed that under conditions of purely centrifugal support 
(neglecting possible electromagnetic corrections of the kind evaluated more re- 
cently [p4|) the condition (33S) for equilibrium, namely the requirement of a 
transcharacteristic rotation speed v = c E , is such that the ring energy is min- 
imised with respect to perturbations preserving the circular symmetry. However 
a more recent investigation of non-axisymmetric perturbations has shown that 
although there are no unstable modes for states of subsonic rotation [Q (as 
exemplified by a cowboy's lassoe loop) with v = c E < a,, instability can nev- 
ertheless occur for rotation in the supersonic regime [f79| , 180 that (contrary 
to what was implicitly assumed in earlier work [|68| , |69| , [ [70| , jn|, 72 



i ffi] using the subsonic type of model given by a linear equation of state for 
which the sum U + T is constant) has been shown by Peter ]7^| , |jT6| , J77J to 
be relevant in the kind of cosmic vortex defects that have been considered so 
far. Although the first category of string loop equilibrium states to have been 
studied systematically has been t hat of circular ring configurations |^3| , it has 
been made clear by recent work [ 1 04 1 , |78| that, as explained above, arbitrary 
non circular equilibrium states are also possible. The stability of such more gen- 
eral equilibrium states has not yet been investigated. While it seems plausible 
that some kinds of "vorton" relic loops may be destroyed by the recently dis- 
covered classical instability mechanism Q , Q , and also the kind of quantum 
tunnelling instability mechanism considered by Davis ]130| ], it docs not seem 
likely that such mechanisms could be so consistently efficient as to prevent the 
long term survival of a lot of other "vorton" equilibrium states. 
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